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ABSTRACT 


Soot  emissions  from  gas-turbine  engines  are  a  concern  for  the  military  for  a  variety  of  reasons 
including  environmental  and  human  health.  A  better  understanding  of  the  process  of  soot 
formation  and  oxidation  can  lead  to  different  strategies  to  mitigate  the  release  of  soot  into  the 
atmosphere.  Studies  have  concentrated  on  soot  formation;  however,  the  oxidation  process  has  not 
been  as  well  characterized,  particularly  the  reactions  that  involve  the  presence  of  O2  which  apply 
to  gas-turbine  engines  which  are  normally  operated  under  excess  O2  concentration.  Soot  oxidizes 
by  reactions  with  molecular  oxygen  (O2),  oxygen  radicals  (O*),  and  hydroxyl  radicals  (OH*). 
The  reactions  of  O*  and  OH*  with  soot  are  relatively  well  understood,  but  the  reactions  of  soot 
with  O2  are  not  as  well  understood.  Therefore,  this  project  focused  on  the  mechanisms  of  soot 
oxidation  by  O2  for  jet  fuels  used  in  military  engines. 

To  separate  the  soot  oxidation  mechanisms  from  the  formation  steps,  a  two-stage  burner  was 
used.  The  two-stage  system  consists  of  an  initial  premixed  burner  where  soot  was  generated 
under  a  variety  of  conditions  for  ethylene/air,  JP-8  surrogate  (n-dodecane/m-xylene)/air,  m- 
xylene/air,  and  n-dodecane/air  premixed  flames.  Downstream,  the  soot-laden  combustion  gases 
were  passed  through  a  secondary  flat-flame  burner  where  soot  was  burned  out  under  fuel  lean  or 
slightly  fuel  rich  conditions.  The  process  of  soot  oxidation  in  the  secondary  burner  was  followed 
by  the  evolution  of  particle  size  distribution  (PSD),  flame  temperature,  gas-phase  composition, 
soot  surface  area,  and  soot  morphology  and  nanostructure  as  a  function  of  the  height  above  the 
second  burner  (HAB).  Measurements  of  soot  size  distribution  and  number  concentration  as  a 
function  of  the  HAB  under  fuel  lean  (Ooverall  =  0.8)  and  slightly  rich  (^overall  =  1.14) 
conditions  showed  particle  fragmentation,  evidenced  by  the  decrease  in  particle  mean  diameter 
and  a  significant  increase  in  number  concentration  in  the  region  where  O2  concentration 
decreased.  Higher  in  the  burner,  soot  was  burned  as  a  result  of  the  increase  in  OH*  concentration 
which  produced  higher  soot  oxidation  rates.  Analysis  of  the  results,  in  terms  of  the  changes  in 
soot  surface  area,  soot  morphology  and  nanostructure  and  the  effectiveness  factor  suggested  soot 
fragmentation  at  a  low  burnout  percentage  with  internal  burning,  which  in  turn,  caused  both,  the 
breakup  of  the  bridges  cementing  primary  particles  and  the  rupture  of  the  primary  particles. 

Experimental  information  from  PSDs,  temperature,  gas-phase  composition  was  used  to  develop 
an  oxidation  kinetic  expression  that  accounts  for  the  effects  of  temperature,  O2,  and  OH*  on  the 
rate  of  oxidation  of  soots  of  different  structure.  The  kinetic  expression  was  able  to  reproduce  the 
major  features  of  the  experimental  results  relative  to  the  other  predictions  using  current  models 
proposed  in  the  literature.  Examination  of  this  additional  aspect  of  oxidation  kinetics  via  O2  will 
enhance  fundamental  combustion  science  and  over  the  long  term  will  enable  military  gas-turbine 
engines  to  continue  to  meet  their  performance  and  operating  requirements  with  reduced  PM 
emissions. 


OBJECTIVES 


SERDP  Relevance 

This  project  responded  to  SERDP’s  Statement  of  Need  (SON)  WPSON-07-01,  Combustion 
Science  to  Reduce  PM  Emissions  for  Military  Platforms.  The  objective  of  this  SON  was  “to 
advance  the  fundamental  combustion  science  relevant  to  non  volatile  particulate  matter  (PM) 
formation  in  military  gas  turbine  engines.”  The  SON  stated  that  modeling  and  experimental 
studies  which  advance  the  capability  to  predict  PM  formation  are  of  interest.  Soot  is  one  of  the 
major  components  of  PM  from  any  combustion  system.  The  formation  of  soot  in  a  gas-turbine 
(or  other)  engine  is  a  complex  process,  involving  soot  inception,  typically  from  acetylene  or 
aromatic  hydrocarbons  such  as  benzene;  soot  growth;  and,  finally  soot  oxidation.  The  project 
focused  on  the  mechanisms  of  soot  oxidation  by  oxygen  for  jet  fuels  used  in  military  engines. 
The  reactions  with  oxygen  are  particularly  important  for  gas  turbines,  which  operate  at  high 
excess-air  rates.  The  examination  of  the  role  of  fuel  components  and  internal  surface  on  soot 
oxidation  rates  was  emphasized  in  the  proposed  study  since  the  allowance  for  these  effects  on  the 
oxidation  rate  are  not  currently  well  characterized.  Examination  of  this  additional  aspect  of 
oxidation  kinetics  enhanced  fundamental  combustion  science  and  over  the  long  term  will  enable 
military  gas-turbine  engines  to  continue  to  meet  their  performance  and  operating  requirements 
with  reduced  PM  emissions. 

Technical  Objectives 

Gas-turbine  engines  are  potentially  a  significant  source  of  PM  for  the  military,  and  PM  emissions 
are  undergoing  increasing  scrutiny  for  their  health  and  regulatory  impacts.  Soot  oxidizes  by 
reactions  with  molecular  oxygen  (O2),  oxygen  radicals  (O*),  and  hydroxyl  radicals  (OH*).  The 
reactions  of  O*  and  OH*  with  soot  are  relatively  well  understood,  but  the  reactions  of  soot  with 
O2  are  not  as  well  studied,  particularly  under  the  conditions  applicable  to  gas-turbine  engines. 
The  proposed  project  focused  on  the  mechanisms  of  soot  oxidation  by  O2  for  jet  fuels  used  in 
military  engines. 

One  method  of  controlling  soot  emissions  is  the  enhancement  of  soot  oxidation.  The  technical 
innovation  of  this  proposal  was  its  focus  on  the  rate  of  soot  oxidation  by  O2.  The  investigators 
hypothesize  that  soot/02  oxidation  kinetics  are  critical  to  the  understanding  of  soot  oxidation  in  a 
gas-turbine  engine.  With  the  high  excess  air  and  the  short  residence  times  in  gas-turbine  engines, 
the  soots  may  have  different  reactivities  than  the  soots  that  are  formed  in  flames  with  higher 
residence  times. 

This  project  had  three  objectives: 

•  To  determine  the  effect  of  the  structure  of  soot,  as  influenced  by  the  fuel  composition  on 
the  rate  of  soot  oxidation  by  O2. 

•  To  quantify  the  role  of  internal  surface  area  on  the  soot  reactivity; 

•  To  develop  power-law  kinetic  correlations  for  soot/02  oxidation  as  a  function  of 
temperature,  O2,  and  time  for  soots  of  different  structures  and  porosity. 
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BACKGROUND 


PM  emissions  from  aircraft  comprise  between  2-8%  of  non-road  mobile  source  emissions  in  the 
United  States,  with  military  aircraft  emissions  responsible  for  13%  of  those  emissions  (EPA 
1991,  1997,  2006;  Blevins  2003).  However,  aircraft  emissions  can  be  one  of  the  major  emission 
sources  on  a  military  base  (Radian,  1998).  Typically,  PM  from  aircraft  is  smaller  than  0.1  pm  in 
diameter  (Champagne,  1971;  Hagen  et  ah,  1998;  Petzold  and  Schroder,  1998;  Pueschel  et  ah, 
1998;  Rogers  et  ah  2005).  Small  particles,  such  as  those  found  in  aircraft  exhaust,  are  thought  to 
deposit  more  efficiently  in  the  lungs  than  larger  particles  and  are  potentially  more  damaging 
(Lighty  et  ah,  2000).  Further,  studies  have  linked  PM  with  an  aerodynamic  diameter  of  2.5  pm  or 
less  (PM2.5)  to  adverse  health  effects  (Dockery  et  ah  1993;  Peters  et  ah  2001;  Tolbert  et  ah 
2000).  Furthermore,  soot  can  be  one  of  the  major  constituents  of  aircraft  emissions  (Rogers  et  ah 
2005),  and  the  soot  content  of  PM  has  been  linked  to  dysrhythmia  (Tolbert  et  ah  2000).  In 
addition  to  health  effects,  aircraft  PM  emissions  are  a  concern  for  climate  forcing  and  cloud 
formation  (Luther  et  ah  1979;  Penner  et  ah  1999). 

Soot  emission  is  the  balance  of  soot  formation  and  oxidation.  There  are  extensive  ongoing 
studies  on  the  difficult  problem  of  soot  formation.  The  proposed  project  focused  on  the  important 
complementary  process  of  soot  oxidation  for  jet  fuels  used  in  military  engines.  In  particular,  the 
effect  of  soot  structure,  both  the  ordering  of  the  carbon  in  the  soot  as  influenced  by  thermal 
annealing  and  the  internal  surface  area,  on  the  rate  of  soot  oxidation  was  studied.  The  reported 
soot  oxidation  rates  by  O2  vary  by  orders  of  magnitude,  probably  because  these  variables 
differed  between  the  various  studies.  It  is  well  known  (e.g.  Shim  and  Hurt  [2000])  that  reactivity 
decreases  with  the  ordering  of  carbons  and  increases  with  increased  internal  surface  area; 
however,  there  have  been  no  systematic  studies  of  these  effects  under  the  temperatures  and  times 
of  interest  to  soot  oxidation  in  gas  turbine  engines. 

Studies  have  shown  that  soot  oxidizes  by  reactions  with  molecular  oxygen  (O2),  oxygen  radical 
(O*),  and  hydroxyl  radical  (OH*).  Fenimore  and  Jones  (1967)  investigated  the  rate  of  soot 
oxidation  under  low  oxygen-gas  partial  pressures  and  temperatures  from  1530  to  1890K  and 
were  the  first  to  report  the  importance  of  OH*  as  a  major  oxidation  reactant.  Their  observations 
suggested  that  about  10%  of  the  soot-OH*  collisions  resulted  in  carbon  removal.  Neoh  et  al. 
(1981)  also  demonstrated  that  OH*  was  the  dominant  oxidant  in  near  stoichiometric  and  rich 
conditions.  The  rate  of  OH*  oxidation  was  described  by  a  collision  efficiency  TOH,  representing 

the  fraction  of  collisions  of  OH*  with  soot  particles  that  resulted  in  the  removal  of  a  carbon 
atom,  times  the  collision  frequency,  i.e. 


W  =  X  £102  t3  FJJtt  Vr  (1) 

where  W  is  the  rate  in  kg/m2*s,  and  P  is  partial  pressure  in  atmospheres.  Neoh  found  the 
collision  efficiency  to  be  0.13  and  relatively  independent  of  equivalence  ratio  and  height  within 
the  flame.  The  literature  suggests  that  OH*  oxidation  is  the  predominant  reaction  mechanism, 
particularly  under  fuel-rich  conditions.  In  addition,  some  oxidation  occurs  via  collisions  with 
O*.  For  example,  Vierbaum  and  Roth  (2003)  found  a  collision  efficiency  of  10%  for  O*.  Table  1 
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show  a  collection  of  collision  efficiencies  for  various  investigators  and  illustrates  that  the 
collision  efficiency  is  relatively  constant  under  various  conditions. 

Table  1.  OH*  collision  efficiencies  for  various  studies. 


Investigators 

Experimental  Conditions 

OH*  Collision  Efficiency 

Bradley  et  al.  1984 

Premixed  flat  flame,  graphite 
carbon 

0.28 

Dai  et  al.  1997 

Laminar  diffusion  flames 

0.14 

El-Leathy  et  al.  2002 

Unsaturated-hydrocarbon/air 
laminar  diffusion  flames 

0.10 

El-Leathy  et  al.  2004 

Acetylene/argon-fueled  high- 
temperature  laminar  diffusion 
flames 

0.13 

Faeth  et  al.  2003 

Acetylene  diffusion  flames  at 
sub-atmospheric  and 
atmospheric  pressures 

0.10 

Fenimore  and  Jones  1967 

Two-stage  burner 

0.10 

Garo  et  al.  1990 

Methane  laminar  diffusion  flame 

0.10 

Kim  et  al.  2004 

Acetylene  laminar  diffusion 
flames  at  pressures  of  0.1 -1.0 
atm 

0.13 

Neoh  1981 

Two-stage  burner 

0.13 

Page  1978 

Two-stage  burner 

0.25 

Puri  et  al.  1994 

Methane/butane  flame 

0.14 

Roth  et  al.  1990 

Diaphragm-type  shock  tube 

0.13-0.31 

Xu  et  al.  2003 

Hydrocarbon/air  diffusion 
flames  at  atmospheric  pressure 

0.14 

The  rate  of  soot  reaction  with  O2  is  not  well  quantified.  Therefore,  the  focus  of  this  work  on 
oxidation  by  O2  was  justified  because  the  rates  of  oxidation  by  OH*  and  O  are  relatively  well 
understood  and  because  the  ratio  of  02/(0H*  +  O*)  is  high  under  the  very  fuel  lean  conditions 
found  in  the  burnout  zone  of  gas  turbines.  The  reaction  rate  of  pyrolytic  graphite  with  O2 
developed  by  Nagle-Strickland-Constable  (NSC)  (1962),  with  parameters  fitted  to  data  obtained 
over  a  temperature  range  of  1273  to  2273K,  is  widely  used  for  soot  oxidation  by  O2.  The  NSC 
relation  was  developed  on  the  assumption  that  the  rate  of  reaction  of  the  graphite  that  they 
studied  was  governed  by  reactions  of  oxygen  with  two  sites:  a  low  reactivity  site  A,  covering  a 
fraction,  x,  of  the  surface,  has  a  rate  governed  by  a  Langmuir-Hinshelwood  expression,  and  a 
high  reactivity  site  B  that  covers  the  remaining  fraction  (1-x)  of  the  surface.  Sites  B  are 
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converted  to  sites  A  in  an  activated  reaction  that  is  weak  function  of  oxygen  concentration.  The 
expression  developed  by  this  model  is: 


w_ 

12 


k  ,  P, 


A1  O, 


1  +  k7P, 


z 1  o 


+  (1  -  x)kBP0 


2  J 


(2) 


where  the  first  term  is  the  contribution  to  the  oxidation  rate  of  the  A  sites  and  the  second  term 
that  of  the  B  sites.  The  fraction  x  of  the  A  sites  are  given  by: 


x  = 


1  + 


kT 

^  a  F’o, 


(3) 


The  parameters  were  fitted  to  the  experimental  data  by  NSC.  These  are  given  in  Table  2  to  show 
the  relatively  high  pre-exponential  constant  and  activation  energy  for  site  A,  which  leads  to  its 
becoming  the  dominant  rate  at  high  temperatures,  and  the  high  activation  energy  for  the 
transformation  of  sites  B  to  A  that  might  be  expected  for  intramolecular  rearrangements. 


Table  2.  Empirical  parameters  for  Nagle  and  Strickland-Constable  Model,  k  =  Ae~E/RT. 


Rate  Constant 

E,  kcal/mole 

A 

Units  for  A 

kA 

30 

20 

g-atom  cm"2  sec"1  atm"1 

kB 

15.2 

4.46x1  O'3 

g-atom  cm"2  sec"1  atm"1 

kj 

97 

1.51xl0'5 

g-atom  cm"2sec"' 

kz 

21.3 

21.3 

atm"1 

Although  the  NSC  relationship  has  been  widely  used  for  soot  oxidation,  it  represents  a  limiting 
low  rate  since  it  applies  to  nonporous  graphitized  carbons.  Graphitic  materials  have  two  types  of 
sites,  the  reactive  edge  carbons  and  the  less  reactive  basal  planes.  Although  NSC  is  an  empirical 
equation,  it  is  consistent  with  such  a  physical  model  with  sites  ‘B’  corresponding  to  the  edge 
sites  and  sites  ‘A’  to  the  basal  plane.  It  is  also  known  that  graphites  can  be  ordered  by  thermal 
annealing,  which  would  correspond  to  the  transformation  of  B  to  A  sites. 

Soots  have  a  wide  range  of  structures  depending  upon  their  temperature  histories.  For  example, 
Du  et  al.  (1990)  have  shown  the  distribution  of  the  sites  through  the  use  of  temperature- 
programmed  desorption.  The  results  shown  in  Figure  1  indicate  that  the  soot  desorption  has  a 
continuous  distribution  of  energies,  rather  than  a  bimodal  distribution.  Hurt  and  Haynes  (2005) 
have  shown  that  such  a  distribution  of  activation  energies  is  common  to  many  carbons  and  that 
they  lead  correlations  for  oxidation  that  show  a  dominant  activation  energy  and  a  power  law 
dependence  on  oxygen. 


5 


10 


c-  8 

a 
o 
x 

\  6 
'o 

E 

w  4 
O 

o 

2 

X  ^ 

S' 

~  o 
-2  . 

20  40  60  80  100 

E  (kcal/mol) 

Figure  1.  Distribution  of  activation  energies  for  oxygen  adsorbed  on  soot  determined  by 
temperature  programmed  desorption  (Du  et  al.,  1990). 

The  structures  of  carbon  vary  from  completely  amorphous  to  complete  ordering  in  planes  in 
graphite;  the  rate  expressions  for  oxidation  will  vary  with  the  extent  of  ordering,  with  the  more 
disordered  structures  having  higher  rate  coefficients.  The  ordering  of  carbon  increases  with  heat 
treatment  and  also  with  extent  of  reaction.  The  ordering  can  be  determined  by  use  of  high- 
resolution  transmission  electron  microscopy  (HR-TEM)  and  image  analysis  as  shown  in  Figure 
2.  The  figure  shows  how  the  ordering  of  a  synthetic  char  (Spherocarb)  increases  with  extent  of 
burnout. 

Another  complicating  factor  is  that  soots  bum  internally  and  develop  high  surface  areas,  as 
shown  in  Figure  3,  which  shows  a  near  tenfold  increase  in  soot  Brunauer-Emmett-Teller  (BET) 
surface  area  with  increases  in  the  extent  of  burnout,  measured  in  a  thermo  gravimetric  analyzer 
(TGA).  The  internal  burning  of  soot  is  important  for  two  reasons:  the  increased  area  leads  to 
higher  reaction  rates,  and  the  increased  porosity  resulting  from  internal  burning  can  lead  to 
particle  fragmentation.  The  variability  of  the  oxidation  rate  of  carbon  with  O2  can  be  appreciated 
from  the  wide  range  of  rates  given  by  the  different  correlations  of  soot  oxidation  summarized  in 
Figure  4,  for  an  assumed  surface  area  of  120  m2/gram  and  an  oxygen  partial  pressure  of  0.1 
atmospheres.  Table  3  provides  selected  information  on  the  experimental  methods  and  carbons 
used  in  the  studies  leading  to  the  rates  summarized  in  Figure  4.  Figure  4  shows  that  the  carbon 
reactivity  varies  by  as  much  as  three  orders  of  magnitude.  In  general,  various  investigators  have 
found  lower  rates  when  using  more  ordered  carbons  for  their  studies.  Clearly  rates  much  higher 
than  those  predicted  by  the  NSC  correlation  have  been  reported.  The  investigators  hypothesize 
that  these  differences  are  a  result  of  different  soot  structures  and  specific  surface  areas  of  the 
soots  or  carbons  used  in  the  different  studies. 


□  Tf=1223K,  11%  conversion 
A  Tj=1223K,  32%  conversion 

O  Tj=1173K.  56%  conversion 
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Figure  2.  Illustration  of  how  the  ordering  of  carbons  varies  with  combustion  conditions,  in  this 
case  increasing  with  carbon  conversion  (Kandas,  1997). 


Figure  3.  Internal  surface  areas  of  two  soots  (ST1  ethylene  soot,  ST2  ethylene  soot  with  2%  Ca 

catalyst)  as  a  function  of  fractional  oxidation. 


Figure  4.  Selected  rates  of  oxidation  of  carbon  by  O2  reported  in  the  literature. 
Table  3.  Studies  used  to  generate  rates  in  Figure  4. 


Number 

Investigators 

Carbons  &  experimental  conditions 

1 

Nagle  and  Strickland-Constable 

High  velocity  oxidant  streaming  over  heated 
graphitized  carbon. 

2 

Park  and  Appleton 

Channel  black  and  furnace  combustion  black,  in 
a  shock  tube  over  1700-3000  K  and  0.05-0.13 
atm  of  oxygen 

3 

Fenimore  and  Jones 

Ethylene-flame-produced  soot,  in  a  two-stage 
burner  over  1530-1800  K  and  0.0001-0.3  atm  of 
oxygen 

4 

Leung  et  al. 

Ethylene  and  propane  laminar  non-premixed 
flames 

5 

Smith 

Porous  and  non-porous  carbons  over  580-2200 
K 

Diesel  soot 

6 

Stanmore  et  al. 

7 

Neeft 

Flame  soot  and  diesel  soot 

8 

Higgens 

Diesel  soot 

9 

Higgens 

Ethylene  diffusion  Flame  soot  generated  in  a 
Santoro-type  burner 

10 

Lee 

Hydrocarbon  laminar  diffusion  flame 

In  most  combustion  carried  out  under  fuel-rich  or  slightly  lean  conditions,  oxidation  by  OH* 
dominates  over  that  by  O2.  Therefore,  few  have  questioned  the  widespread  use  of  the  NSC  rate 
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for  O2  since  it  was  a  minor  contributor  to  soot  oxidation.  But  the  NSC  rate  with  its  two  reactive 
sites  is  best  used  for  graphitic-like  structures  that  have  predominantly  two  carbon  sites  and 
intrinsically  lower  reactivities.  Studies  on  chars  with  disordered  structures,  by  contrast,  have 
shown  that  carbon  reactivity  can  be  modeled  by  sites  with  a  distribution  of  adsorption  and 
desorption  energies  (Haynes  and  Newberry  2000).  This  approach  leads  to  a  rate  expression  with 
a  fractional  order  for  gaseous  oxygen  (Hurt  and  Haynes  2005),  rationalizing  the  power-law 
correlations,  which  have  been  used  for  decades.  It  is  the  objective  of  the  proposed  study  to 
investigate  the  oxidation  of  soot  produced  under  conditions  relevant  to  gas-turbine  combustors 
and  to  develop  the  appropriate  kinetic  expressions. 

MATERIALS  and  METHODS 

The  study,  designed  to  take  advantage  of  the  developments  in  sampling,  sizing,  and  structural 
characterization  of  soot,  investigated  the  kinetics  of  soot  oxidation  with  O2  and  OH*  and 
examined  the  role  of  soot  structure,  as  influenced  by  fuel,  soot  surface  area  and  particle  size 
distribution. 

Burner  System 

The  University  of  Utah’s  two-stage  burner  (Figure  5),  designed  after  the  apparatus  of  Neoh 
(1981),  isolates  the  soot  oxidation  mechanisms  from  the  formation  steps  (Neoh,  1981;  Merrill 
2005;  Lighty  et  al.,  2008,  Echavarria  et  al.,  2011).  In  this  system,  soot  was  generated  in  a  fuel- 
rich  premixed  flame  (bottom  burner),  which  serves  as  the  first  stage.  Soot  was  then  oxidized 
under  lean  conditions  and  slightly  rich  conditions  in  the  top  burner.  Soot  was  derived  from 
different  sources  to  account  for  the  effect  of  the  parent  fuel  on  soot  oxidation.  In  the  first 
experiments  soot  was  produced  in  an  ethylene/air  flame  (Oi  =2.5).  For  the  jet  fuels  or  surrogate 
mixtures  soot  was  formed  in  a  surrogate/air  flame  (<E>i  =  2.11).  The  surrogate  consists  of  a 
mixture  m-xylene/n-dodecane  with  a  volume  ratio  of  10/90  (referred  to  as  “surrogate”)  which 
approximates  the  range  of  aromatic  concentrations  in  JP-8  (Edwards,  2005).  For  comparison, 
soot  was  also  generated  using  m-xylene/air  (Oi  =  1.7)  and  n-dodecane/air  flames  (d>i  =2.15)  that 
represented  the  pure  components  of  the  surrogate  mixture.  The  soot-laden  combustion  gases 
were  passed  through  a  secondary  flat-flame  premixed  burner  where  soot  was  burned  out  under 
fuel  lean  (Ooverall  =  0.8)  and  slightly  fuel  rich  (Ooverall  =  1.14)  conditions.  The  overall 
equivalence  ratio  is  defined  by  the  following  equation: 

9* overall  =  [Fuel/(Airi+Air2)]actual/[Fuel/Air]stoichiometric,  (4) 

and  the  stoichiometry  of  the  secondary  burner  was  controlled  by  the  injection  of  air  and/or  H2 
(see  Table  4)  through  the  secondary  air  port.  H2  helped  to  obtain  more  stable  flames  in  the  top 
burner  for  the  conditions  where  soot  was  generated  in  m-xylene/air  and  n-dodecane/air  flames. 
Air  and  gaseous  fuels  were  fed  into  the  two-stage  burner  using  5850  E  mass  flow  controllers  and 
liquid  fuels  (surrogate)  were  vaporized  and  mixed  with  preheated  air  prior  to  entering  the  first 
burner.  Temperature  along  the  feeding  line  and  below  the  tube  bundle  in  the  first  burner  was 
controlled  at  250  °C  to  avoid  fuel  recondensation.  An  initial  vaporization  system  consisted  of  a 
bubbler  which  was  later  changed  for  a  commercial  (Mesoscopic  Devices 
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Air 


Table  4.  Summary  of  Experimental  Conditions 


Bottom  Burner 

Top  Burner  (secondary  gas  inlets) 

System  O  overall  Fuel,  g/min  Air,  g/min 

Air,  g/min  H?,  g/min 

,  ...  0.8  0.734  4.322 

Ethylene/  Air 

J  1.14  0.734  4.322 

9.220  0 

5.135  0 

c,  .  ...  *  0.8  0.835  5.844 

Surrogate  /Air* 

1.14  0.835  5.844 

9.570  0 

4.960  0 

v  ,  ...  0.87  0.728  5.797 

m-Xylene/Air 

1.14  0.728  5.797 

7.408  0.0466 

4.280  0.0466 

,  ...  0.8  0.843  5.852 

n-Dodecane/Air 

1.14  0.843  5.852 

10.945  0.0254 

5.941  0.0254 

*  Surrogate  mixture  consist  of  10  %  volm-Xylene  and  90  %vol.  n-Dodecane 


Inc.)  vaporizer  coupled  to  a  syringe  pump  and  temperature  control  system  (see  Figure  6).  Both 
systems  provided  a  uniform  fuel  flow  to  the  burner.  However,  the  system  with  the  commercial 
vaporizer  allowed  a  wider  range  of  flow  conditions  and  ensured  that  the  vaporization  process  did 
not  alter  the  fuel  through  thermal  chemical  reaction  or  preferential  vaporization.  To  verify  this, 
the  vaporized  fuel  was  captured  in  a  cold  solvent  (dichloromethane)  trap  and  analyzed  using  gas 
chromatography.  The  concentration  of  liquid  samples  before  vaporization  and  after  vaporization- 
condensation  was  found  to  be  the  same  (see  Figure  7).  Experimental  measurements  of 
temperature,  particle  size  distribution,  gas-phase  composition,  soot  surface  area,  and  soot 
morphology  and  nanostructure  were  mainly  performed  within  the  top  burner. 
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Figure  6.  Schematic  representation  of  the  bottom  burner  showing  injection  system  for  liquid 

fuels. 
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Figure  7.  Chromatograms  of  the  surrogate  mixture  before  vaporization  and  after  vaporization- 

condensation. 


Flame  Temperature  Measurements 

Temperature  profiles  along  the  centerline  of  the  top  burner  were  measured  using  a  0.008” 
diameter  uncoated  Type-B  thermocouples.  Comparisons  to  uncoated  versus  coated  yielded  a 
difference  in  temperature  of  approximately  40°C  which  is  below  the  typical  uncertainty  (±50- 
100°C)  for  this  thermocouple  (McEnally  et  al.,  1997;  Shaddix,  1999).  The  thermocouple  was 
inserted  into  the  flame  using  a  fast  insertion  mechanism  (Figure  8).  In  this  system  the  transient 
response  of  the  thermocouple  was  recorded  for  25  seconds  at  a  sample  rate  of  50  samples  per 
second.  According  to  McEnally,  et  al.  (1997),  when  a  thermocouple  is  inserted  into  a  sooty 
flame,  particles  will  deposit  on  its  junction  due  to  the  thermal  gradient  between  the  gas  and  the 
thermocouple.  Thus  a  correlation  of  the  transient  response  and  the  soot  deposition  rate  could  be 
used  to  estimate  the  local  gas  temperature  and  the  local  soot  volume  fraction  without  having  to 
make  assumptions  about  particle  size  distribution,  particle  optical  properties  or  gas  uniformity. 


Figure  8.  Fast  response  thermocouple  setup. 


Figure  9  presents  typical  profiles  of  the  transient  response  of  the  thermocouple  for  both 
nonsooting  and  sooting  conditions  measured  in  the  bottom  burner  for  ethylene/air  flames.  Under 
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nonsooting  conditions,  the  temperature  of  the  bead  increased  rapidly  and  leveled  off  after 
reaching  a  maximum  (see  Figure  9a).  By  contrast,  when  soot  is  in  the  flame,  it  was  driven  to  the 
bead  surface  because  of  the  thermophoretic  gradient  between  the  flame  environment  and  the  cold 
bead  surface.  As  result,  temperature  dropped  continuously  as  soot  deposited  on  the  thermocouple 
bead,  increasing  both  its  diameter  and  emissivity.  Once  the  bead  was  sufficiently  coated,  the 
emissivity  of  the  bead  was  equal  to  the  emissivity  of  the  soot,  and  the  temperature  kept 
decreasing  but  at  a  different  rate.  In  the  top  burner,  temperature  profiles  showed  similar  behavior 
relative  to  Figure  9a. 
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Figure  9.  Transient  response  of  the  thermocouple  under: 
(a)  Nonsooting  conditions  and  (b)  Sooting  conditions. 
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The  temperature  of  the  gases  ( Tg )  can  be  inferred  by  a  quasi-steady  energy  balance  at  the 
junction  as  follows: 


2  d, 


(n-n) 


(5) 


where  7}  is  the  junction  temperature,  gj  is  the  junction/bead  emissivity,  cr  is  the  Stefan-Boltzmann 
constant,  Niij  is  the  junction  Nusselt  number,  dj  is  the  junction  diameter,  and  kgo  is  the  gas 
thermal  conductivity.  7}  can  be  estimated  from  the  response  curve  of  the  thermocouple  (Figure  9) 
by  extrapolating  the  temperature  at  t  =  0  seconds.  For  this  case,  where  we  used  B-type 
thermocouples,  the  emissivity  can  be  estimated  as  a  function  of  the  temperature  as  follows: 

s  =  0.1083  ln(Ty )  -  0.5644  (6) 

The  Nusselt  number  has  been  evaluated  for  similar  systems  in  the  range  from  2.26  to  2.35.  kgo  is 
assumed  constant,  which  is  reasonable  for  combustion  gases  at  high  temperature.  Radiation 
correction  performed  for  the  data  in  Figure  9  produced  Tg  of  1648  K  and  1580  K  for  the 
nonsooting  and  sooting  conditions  respectively. 

Particle  Size  and  Concentration 

Particle  size  distributions  and  soot  concentration  were  measured  with  a  scanning  mobility 
particle  sizer  (SMPS).  This  system  consisted  of  a  TSI  3080  classifier  with  nano  and  long 
differential  mobility  analyzers  (DMA)  and  a  3025  ultrafme  condensation  particle  counter 
(UCPC).  Nano  and  long-DMAs  and  different  impactor  sizes  (0.0457  and  0.071  cm)  were  used  to 
cover  PSDs  in  the  size  range  from  3  to  660  nm  (SMPS  User’s  Manual).  A  general  overview  of 
the  SMPS  and  sampling  system  is  presented  in  Figure  10.  In  this  system,  particles  from  the  flame 
were  drawn  into  a  horizontal  sampling  probe  through  a  small  orifice  (0.32  mm  diameter)  and 
instantaneously  diluted  using  a  N2  stream  at  30  1pm  (STP).  A  small  portion  (1  1pm)  of  the  exit 
stream  was  sent  to  the  SMPS  for  particle  sizing  and  number  concentration  measurements.  The 
large  portion  of  the  exit  was  additionally  diluted  with  air  in  an  eductor  and  vented  through  a 
hood.  The  eductor  not  only  diluted  the  large  stream  but  also  allowed  controlling  the  pressure 
drop  through  the  pinhole. 

In  the  DMA  (Chen  et  al.,  1998;  Maricq  et  ah,  2000,  2004,  2005),  an  electrostatic  classifier 
extracts  a  particular  size  fraction  of  particles  from  the  polydisperse  stream.  The  size-selected 
particle  stream  enters  a  Kr-85  neutralizer  where  the  particles  are  engaged  in  frequent  collisions 
with  bipolar  ions.  Equilibrium  is  quickly  reached  and  the  particles  carry  a  bipolar  charge 
distribution.  This  charged  aerosol  then  travels  to  the  differential  mobility  analyzer.  The  DMA 
consists  of  two  concentric  cylinders:  a  collector  rod  that  is  maintained  at  a  set  negative  voltage 
and  an  electrically  grounded  outer  cylinder.  This  cylinder  orientation  creates  an  electric  field. 
This  electric  field  leads  to  the  attraction  of  positively  charged  particles  to  the  negatively  charged 
collector  rod.  Particles  within  the  user-defined  range  of  electrical  mobility  exit  the  collector  rod 
through  a  small  slit  and  travel  to  an  Ultrafme  Condensation  Particle  Counter  (UCPC)  where  the 
particles  pass  over  a  heated  pool  of  alcohol  and  are  saturated  with  alcohol  vapor.  Next,  the 
aerosol  passes  into  a  condenser  and  is  cooled.  The  alcohol  condenses  onto  the  surface  of  the 
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particles  and  they  reach  a  size  that  is  optically  visible.  At  this  point,  the  particles  can  be  counted. 
The  flow  schematic  for  the  electrostatic  classifier  is  shown  in  Figure  1 1 . 
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Figure  10.  SMPS  and  dilution  sampling  system  indicating  probe  ID,  OD  and  line  lengths. 


Pump 


Figure  11.  Flow  Schematic  for  the  Electrostatic  Classifier  (TSI  Inc.). 

Dilution  During  the  Sampling  of  Nanoparticles  from  Flames 

An  essential  part  in  the  study  of  aerosols  is  the  ability  to  collect  representative  samples  for 
analysis.  These  samples  must  accurately  reflect  the  aerosol  concentration  and  the  particle  size 
distribution.  During  online  analysis  by  the  SMPS  in  a  flame,  sampling  issues  range  from  particle 
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coagulation  and  agglomeration  to  diffusive  losses  in  the  sampling  probe  and  sampling  line.  The 
most  common  method  to  minimize  these  problems  during  sampling  from  flames  is  by  diluting 
the  soot-laden  combustion  gases.  The  use  of  a  cold  dilution  gas  leads  to  an  immediate  quench  of 
particle  growth  chemistry  and  minimizes  the  thermophoretic  deposition  of  soot  in  the  sample  line 
that  occurs  when  high-temperature,  soot-laden  gases  come  in  contact  with  a  cold  surface.  Zhao  et 
al.  (2003)  showed  that  particle  diffusion  losses  and  particle  coagulation  can  be  minimized  by 
systematically  increasing  the  dilution  ratio  to  a  critical  value  where  the  particles  size  distribution 
function  becomes  independent  of  the  dilution  ratio.  Experimentally,  this  critical  dilution  ratio 
was  around  104.  However,  other  studies  have  shown  that  this  value  can  range  from  103  to  104 
depending  on  equivalence  ratio  and  combustion  system  utilized  (Maricq  et  al.,  2001;  Zhao  et  al., 
2003;  Manzello  et  al.,  2007;  Echavarria  et  al.,  2009;  Abid  et  al.,  2009;  D’Anna,  2009). 

In  this  study,  we  used  a  dilution  system,  similar  to  that  of  Zhao  et  al.  (2003)  and  Kasper  et  al. 
(1997)  with  changes  in  probe  size  (OD  =  9.51  mm)  and  pinhole  diameter  (0.32  mm),  yielding 
dilution  ratios  greater  than  104,  which  minimized  wall  losses  and  quenched  reactions  or 
coagulation  that  would  otherwise  occur  in  the  sampling  system.  Although,  this  system  minimized 
diffusion  losses  and  particle  coagulation  during  sampling,  corrections  are  still  necessary  to 
ensure  that  the  final  result  reflects  the  actual  evolution  of  the  particle  size  in  our  flames. 
Corrections  for  penetration  efficiency,  into  the  probe  and  probe  orifice,  and  diffusion  losses 
during  transport  were  mainly  applied  following  the  procedure  presented  by  Minutolo  et  al. 
(2008).  Corrections  due  to  diffusion  losses  in  the  SMPS  were  conducted  using  the  AIM  software 
upgrade  (TSI,  In.,  2006). 

Dilution  Ratio  Calculation 

As  previously  mentioned,  high  dilution  in  the  sampling  point  was  critical  to  quench  and 
minimize  chemical  reactions,  coagulation  and  wall  particle  loses  along  the  sampling  system  (see 
Figure  10).  The  dilution  ratio  (DR)  was  defined  as  the  volumetric  flow  of  the  dilution  gas  (N2,  30 
lpm,  STP)  divided  by  the  sampled  flow  through  the  pinhole  (Qp,  STP),  or 


Qp(STP)  (7) 

Qp  was  corrected  to  standard  temperature  and  pressure  (STP)  conditions  from  the  temperature 
and  pressure  condition  at  the  pinhole  orifice  (~  1000  K  ,  Patm).  With  this  definition,  DR  can  be 
evaluated  as  a  function  of  the  pressure  drop  ( APp )  or  flow  (Qp)  through  the  pinhole  orifice.  Qp 
was  modified  by  adjusting  the  pressure  drop  (API  and  API)  through  the  sampling  probe, 
accomplished  via  an  eductor  (Fox  Valve  Mini-eductors,  Model  611210-060)  installed 
downstream  from  the  probe. 

The  eductor  (Figure  12)  was  used  similarly  to  a  vacuum  pump.  However,  instead  of  using 
electricity,  it  uses  fluid  mechanics  to  create  certain  amount  of  suction  from  the  energy  that  is 
stored  in  a  motive  gas  (air).  The  motive  fluid  was  passed  through  a  precision-machined  nozzle. 
This  led  to  a  high  velocity  jet  that  induced  a  vacuum  that  pulled  the  diluted  sample  through  the 
suction  port.  The  inlet  differential  pressure  (AP3)  of  the  motive  gas  was  adjusted  to  increase  or 
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decrease  the  suction,  which  in  turn,  changed  the  pressure  drop  through  the  probe  and  pinhole 
orifice. 
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Discharge 


Suction  A 

Figure  12.  Schematic  representation  of  the  eductor. 


Zhao  et  al.,  (2003)  showed  that  the  dependence  of  Qp  on  the  pressure  drop  ( APp ),  pinhole 
diameter  (dp)  and  length  (/)  can  be  described  by  the  hydrodynamic  theory  for  a  low  Reynolds 
number  as  follows: 


Qp  =  c 


71  dp* 
Ip 


(8) 


where  c  =  0.0015,  and  the  viscosity  (p)  was  obtained  using  the  premixed  code  in  CHEMKIN  4.1 
for  each  condition,  dp  was  equal  to  0.24  mm  and  the  length  of  the  orifice  was  equal  to  0.65  mm. 
These  values  were  the  result  of  a  series  of  preliminary  tests  using  different  pinhole  orifice  sizes 
in  order  to  minimize  particle  losses  and  coagulation  in  the  dilution  point.  The  pinhole  diameter 
(dp)  was  found  to  be  satisfactory  in  the  range  from  0.2  to  0.32  mm  ID.  Below  this  range, 
significant  particle  losses  and  pinhole  clogging  was  observed,  and  above  0.32  mm  it  was  not 
possible  to  get  enough  dilution  to  minimize  wall  losses  and  coagulation  along  the  sampling  line. 
The  length  of  the  pinhole  orifice  (/)  was  reduced  until  the  PSD  converged  to  a  single  result.  The 
pressure  drop  through  the  pinhole  orifice  (APp)  was  estimated  as  the  average  between  API  and 
AP2.  API  and  AP2  were  measured  experimentally  using  U  manometers,  since  the  pressure  drop 
was  expected  to  be  very  small;  an  oil-like  fluid  with  a  specific  gravity  equal  to  0.856  was 
utilized.  To  use  Equation  (7)  and  (8)  correctly,  it  was  necessary  to  account  for  the  change  in 
pinhole  diameter  due  to  soot  deposition. 

Figure  13  presents  the  change  Qp  through  the  pinhole  as  a  function  of  AP3  and  the  change  in 
dilution  ratio  with  Qp  for  studies  on  soot  oxidation.  Qp  was  observed  to  increase  with  increasing 
the  inlet  pressure  through  the  eductor  (AP3),  which  was  a  direct  result  of  the  increase  in  suction. 
On  the  other  hand,  the  dilution  ratio  was  found  to  be  inversely  proportional  to  Qp  as  expected 
and  as  shown  in  other  studies  (Zhao  et  al.,  2003,  Abid  et  al.,  2008).  Experiments  during  soot 
oxidation  were  mainly  carried  out  in  the  DR  range  from  103  to  104.  The  uncertainty  of  this 
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methodology  to  calculate  the  dilution  ratio  has  been  reported  in  the  literature  in  the  range  from 
20  to  40  %  (Manzello  et  ah,  2007;  Echavarria  et  ah,  2009;  Abid  et  al.,  2009).  The  experimental 
error  found  in  this  study  fell  into  that  range  and  it  became  the  main  source  of  uncertainty  in  the 
measurements  of  number  and  mass  concentration  obtained  from  the  integration  of  the  PSDs. 
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Figure  13.  Dilution  ratio  calibration  (a)  Qp  as  a  function  of  AP3  and  (b)  DR  as  a  function  of  Qp. 


Corrections  to  PSDs 

Corrections  to  PSDs  for  penetration  efficiency  into  the  probe  and  probe  orifice,  diffusion  losses 
during  transport  and  diffusion  losses  in  the  SMPS  were  necessary  to  ensure  that  the  final  result 
reflected  the  actual  evolution  of  the  particle  size  in  our  flames.  Besides  coagulation,  particle 
diffusion  losses  along  the  sampling  line  and  in  the  SMPS  occur  when  the  particles  collide  with  a 
surface  and  remained  together  due  to  Vander  Wal  forces,  electrostatic  forces,  and  surface  tension 
(Hinds,  1982).  Another  mechanism  that  can  favor  the  diffusion  of  particles  to  the  walls  of  the 
probe  is  the  thermophoresis.  However,  it  becomes  negligible  in  this  case  since  the  temperature 
gradient  inside  the  dilution  probe  was  very  small.  PSDs  are  usually  presented  in  the  form  of 
dN/dlogDp  versus  Dp.  Even  though  the  SMPS  can  provide  results  in  this  form,  it  does  not 
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account  for  particle  losses  and  coagulation  along  the  sampling  line;  therefore,  it  is  necessary  and 
recommended  to  recalculate  dN/dlogDp  from  the  raw  counts  provided  by  the  SMPS  as  follows: 

dN,  _  Ci*DR 

d  log  Dp  Q*  tSMPS  *  tj  *  d  log  Dp)  (9) 

where  Ci  are  the  raw  counts  provided  by  the  SMPS,  DR  is  the  dilution  ratio,  and  Q  and  (,smps  are 
the  sample  flow  and  sample  time  through  the  SMPS,  respectively.  The  product  Q*tsMPS 
represents  the  total  volume  of  air  sampled,  r/  is  the  total  sampling  efficiency  which  is  given  by 
the  product  of  the  SMPS  efficiency,  efficiency  into  the  probe,  and  efficiency  through  the 
sampling  line. 


V  =  V SMPS  *V probe  *V line 


(10) 


The  SMPS  efficiency  (  tjsmps )  can  be  estimated  as  the  product  of  the  efficiency  for  the  ionization 
process  (>/„„),  and  the  diffusion  losses  correction  inside  the  instrument  rjsMPs-diff-  The  ionization 
efficiency  for  a  similar  instrument  has  been  reported  by  Sgro  et  al.  (2008,  2009).  It  was  based  on 
Wiedensohler's  equation,  which  describes  the  fraction  of  charged  particles  with  diameters 
1  <£)/?<  1000  nm.  The  penetration  efficiency  through  the  SMPS  ( rjsMPs-diff )  has  been  estimated 
using  the  diffusion  losses  correction  application  of  the  SMPS  software  upgrade  (TSI,  Inc,  2006). 
This  application  includes  the  penetration  efficiencies  along  the  SMPS  system. 


The  penetration  efficiency  into  the  sampling  probe  ( rjprobe )  accounts  for  the  sampling  efficiency 
into  the  line  h  where  the  flow  is  turbulent  and  in  the  sampling  orifice  where  the  Reynolds 
number  may  be  approximated  as  laminar  or  creeping  flow.  The  turbulent  sampling  efficiency 
( r/r )  can  be  estimated  by: 


ijT  =  exp 


' -7t* D* /j  * Sh^ 

Qn2  j 


(11) 


Sh  =  0.01 18*  Re7/8ScI/3 


(12) 


where  D  is  the  diffusion  coefficient  of  particles,  /;  is  the  length  of  the  sampling  probe,  QN2  is  the 
volumetric  flow  rate  of  the  dilution  gas  corrected  for  the  temperature  and  pressure  conditions 
inside  probe,  Sh  is  the  Sherwood  number  and  is  given  by  Equation  (12).  The  diffusion  coefficient 
can  be  estimated  using  the  Stokes-Einstein  equation  for  aerosols  in  terms  of  the  particle  mobility 
as  follows  (Hinds,  1982): 


D  = 


k*T*Cc 
3n*ju*Dp  ■ 


(13) 


where  k  is  the  Boltzmann  constant,  T  is  the  temperature  inside  the  probe  which  was  found  to  be 
approximately  373  K,  //  is  the  dynamic  viscosity,  Dp  is  the  particle  diameter  and  C  is  the 
Cunningham  slip  correction  factor  calculated  by  (Kim  et  al.,  2004): 
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(14) 


C(Kn)  =  \  +  A*  Kn 


A  -  a  +  /3*  exp 


with  parameters  a 
a  Knudsen  number  from  0.5  to  83. 


zL' 

\Kn  ) 

1.165,  =  0.483,  and  y 


(15) 

0.997  fitted  for  particle  sizes  up  to  270  nm  and  for 


To  calculate  the  penetration  efficiency  through  the  pinhole  orifice  ( rjp)  two  different  approaches 
were  used.  First,  flow  through  the  pinhole  orifice  was  assumed  to  be  laminar  and  the  penetration 
under  this  condition  is  given  by  (Minutolo  et  al.,  2008): 


( LP  -  0.0191  *  esp(-S.eS7  if  "{•  3  +  0.097S  *  b«:jC-22,3  if  V )  +  0,03250  *  eaX-B?  *  ■(■  1 

fQT  (£  0,0312, 

&  QF  -  1  -  2,56(7(2/8)  +  1.2  #  £  +Q.177£T(4/8)  for  £«  0.0812  (16) 

n  *D*l 


5  = 


Qp 


(17) 


where  5  is  the  deposition  parameter  and  1  is  the  length  of  the  tube.  The  particle  diffusion 
coefficient  can  be  calculated  using  Equation  (13)  for  the  temperature  of  the  combustion  products 
entering  the  pinhole  orifice  (-1000K).  A  similar  approach  for  penetration  efficiency  in  a  tube 
under  laminar  conditions  can  be  found  in  Hinds  (1982). 

It  is  essential  to  note  that  there  were  several  uncertainties  associated  with  the  method  described 
above  to  correct  for  particles  losses  inside  the  sampling  orifice  (Sgro  et  al.,  2009).  First, 
Reynolds  numbers  were  between  0.3  to  25  depending  on  the  dilution  ratio,  which  suggests  that 
the  flow  should  be  consider  as  creeping  flow  instead  of  laminar.  As  result,  particle  losses  may  be 
higher  or  lower  than  those  estimated  from  laminar  flow  theory  which  assumes  fully  develop 
flow.  Second,  the  temperature  at  the  sampling  orifice  was  extrapolated  rather  than  measured 
experimentally.  It  has  a  direct  effect  on  properties  such  as  the  particle  diffusion  coefficient, 
viscosity,  density  etc.  which  can  impact  the  final  results  considerably.  Even  though  these 
uncertainties  can  affect  the  final  result,  significant  efforts  were  performed  during  the  experiments 
to  minimize  them. 

Another  approximation  to  estimate  particle  losses  in  the  sampling  orifice  has  been  discussed  by 
Sgro  et  al.,  (2009).  In  this  approach,  the  particle  sticking  efficiency  (/(Dp))  is  evaluated  as  a 
function  of  the  particle  diameter  by  the  interaction  potential  between  the  particles  and  the  orifice 
walls.  Additionally,  the  effect  of  particle  rebound  after  particles  collide  with  the  walls  can  be 
accounted  for  assuming  that  the  number  of  collisions  with  the  walls  producing  particle  losses  is 
equal  to  the  product  of  the  number  of  collisions  and  the  particle  sticking  efficiency.  Even  though 
this  second  approach  seems  to  be  more  realistic,  no  significant  difference  was  obtained  either  by 
using  it  or  the  assumption  of  laminar  flow  through  the  pinhole  orifice. 

Figure  14  presents  a  summary  of  a  typical  result  for  rjion,  t/smps,  *1t,  Pane,  ?]P,  and  the  overall 
penetration  efficiency  (77).  It  should  be  taken  into  account  that  the  total  penetration  efficiency  can 
be  considerably  affected  by  different  factors  such  as  the  dilution  ratio  used  to  collect  the  PSD 
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and  the  length  of  the  sampling  line  in  the  laminar  side  (I2+I3).  The  higher  the  dilution  ratio  the 
lower  the  flow  through  the  pinhole,  which  in  turn,  reduced  the  sampling  efficiency  of  the  smaller 
particles  and  increased  the  experimental  error.  In  addition,  shortening  the  sampling  line  reduced 
particle  losses,  particularly  in  the  laminar  side.  Using  the  data  provided  by  Figure  14,  the 
corrected  PSDs  can  be  obtained  from  the  raw  counts  obtained  in  the  SMPS. 


Figure  14.  Penetration  efficiencies 
thorough  the  turbulent  side  (r|T),  laminar 
side  (111™.),  sampling  orifice  (r|P), 
SMPS(t|smps)  and  overall  sampling 
efficiency  ( r\ ). 


Chemical  Analysis 

H2,  O2,  CO  and  CO2  were  isokinetically  sampled  and  analyzed  by  online  gas  chromatography 
(VARIAN,  CP-4900  Micro  GC).  The  sampling  system  (Figure  15)  consist  of  a  water-cooled 
probe  (0.2  cm  ID)  coupled  to  a  filter  ( 1  pm)  trap  followed  by  a  chiller  where  water  is  removed 
from  the  gaseous  stream.  Flow  through  the  system  is  controlled  by  a  rotameter  located  in  the 
inlet  of  the  vacuum  pump.  A  small  sample  of  the  exit  stream  is  sent  to  the  GC  for  gas-phase 
analysis  and  the  other  part  is  diluted  with  N2  and  vented.  The  temperature  inside  the  probe  and  in 
the  filter  trap  was  kept  between  150  -  200  °C  to  avoid  water  condensation. 
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Figure  15.  Sampling  system  to  measure  gas-phase  compounds. 
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Soot  Surface  Area  Measurements 

The  surface  area  of  the  soot  sample  was  measured  gravimetrically  by  adsorption  of  CO2  at  297K. 
The  method  consisted  of  injecting  increasing  concentrations  of  high  purity  CO2  in  helium  (He) 
(0,  0.1,  0.2,  0.35,  0.5  and  1  mole  fractions)  into  the  TGA,  and  recording  the  weight  change  due  to 
CO2  adsorption.  Initially,  after  the  sample  has  been  loaded  into  the  TGA,  the  system  was  heated 
up  to  150  °C  in  He  or  N2  to  remove  adsorbed  water.  When  the  system  cooled  down,  He  was 
injected  into  the  system  at  1.85  1pm  (STP).  Between  40  to  60  minutes  were  allowed  for  the 
system  to  come  to  equilibrium  for  each  point.  Afterwards,  CO2  was  fed  into  the  TGA  at  a  molar 
fraction  of  0. 1  in  He.  Four  additional  He/CCE  gas  mixtures  switches  were  performed  while  the 
total  flow  and  the  system  pressure  (Patm)  were  maintained  constant. 

Buoyancy  and  drag  effects  are  usually  a  concern  during  measurements  of  weight  change  in 
TGAs.  To  correctly  estimate  the  mass  of  CO2  adsorbed  in  the  solid  samples  these  effect  must  be 
accounted  for.  These  effects  are  associated  with  apparent  loss  or  gain  in  weight  depending  on 
changes  of  density  in  the  system.  A  blank,  where  an  experiment  was  run  under  the  same 
experimental  conditions  with  either  an  empty  sample  holder  or  the  holder  with  an  inert  material 
(like  quartz),  was  carried  out  to  account  for  the  weight  changes  due  to  buoyancy  and  drag  effects 
(Campbell,  2005). 


To  test  this  methodology  a  powder  coal  sample  (CA  Corp.)  with  known  BET  surface  area  was 
used.  The  BET  surface  area  for  this  sample  was  determined  using  a  Tristar  II  3020  Surface  Area 
Analyzer  (Micromeritics).  The  value  found  for  this  sample  using  the  surface  area  analyzer  was 
660  ±  10  m2/g.  Next  the  sample  was  run  using  the  methodology  described  above.  Figure  16 
presents  a  complete  CO2  adsorption  experiment.  During  the  two  first  injections  of  CO2  into  the 
system  (Yco2  =  0.1  and  Yco2  =  0.2),  there  were  slight  increases  in  weight  which  were  associated 
with  CO2  adsorption  on  the  sample.  For  the  other  points  (Y<x>2  =  0.35,  0.5  and  1)  adsorption  is 
still  occurring,  however,  weight  dropped  due  to  the  major  contribution  of  buoyancy  effects 
caused  by  the  difference  in  the  densities  of  He  and  CO2.  At  the  end  of  the  experiment,  CO2  was 
again  replaced  by  He,  and  CO2  desorbed  from  the  sample. 


Figure  17  illustrates  a  test  run  under  the  same  experimental  conditions  with  an  empty  sample 
holder.  This  result  was  an  indication  that  all  the  weight  changes  can  be  associated  with  CO2 
adsorption  and  to  drag  and  buoyancy  effects,  which  at  higher  CO2  concentrations  accounted  for 
the  major  weight  changes.  These  measurements  with  and  without  sample  allowed  establishing  a 
relationship  between  the  amount  of  CO2  adsorbed  and  the  CO2  concentration  or  the  partial 
pressure  in  the  system.  Figure  18  presents  adsorption  isotherms  ( AW/Wo  versus  Pco2,  Wo  is  the 
initial  sample  weight)  at  297  K  for  the  commercial  coal  powder,  and  the  soot  samples  collected 
in  the  two- stage  burner  at  0  and  2.5  mm  under  ^overall  =  0.8.  The  methodology  to  extract  the 
results  for  the  soot  samples  was  the  same  to  that  in  the  coal  sample  (CA  corp.)  These  results 
were  further  analyzed  to  obtain  soot  surface  area  using  the  Dubinnin-Polayi  (DP)  method  (Du, 
1990;  Kandas,  1997).  The  DP  method  used  Equations  (18)  and  (19)  to  estimate  the  CO2 
micropore  adsorption  capacity  of  the  samples: 
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Figure  16.  CO2  adsorption  test  on  commercial  coal  powder  (CA  Corp.)  at  297  K. 


Figure  17.  CO2  adsorption  test 
carried  out  on  an  empty  sample 
holder.  Weight  changes  were 
due  to  drag  and  buoyancy 
effects. 
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A  plot  of  log  x  versus  log  (Ps/Pcoi)  yielded  straight  lines  (see  Figure  19)  for  the  coal  sample 
and  for  the  two  soot  sample  collected  in  the  two- stage  burner.  The  monolayer  adsorption 
capacity  was  obtained  from  the  intercept  and  the  specific  micropore  surface  area  was  calculated 
as  follows: 


SA(™2  I  g sample  ) 


xsNocj 


co2 


MC02 


(20) 


No  is  the  Avogadro's  number,  oco2  is  the  cross-sectional  area  of  a  CO2  molecule  (0.17  nm2)  and 
Mco2  is  the  CO2  molecular  weight.  Using  this  approach,  surface  areas  of  560  ±  40  m2/g  ,  282  ± 
35  m2/g  and  494  ±  27  m2/g  were  obtained  for  the  coal  powder  sample,  soot  sample  at  0  mm 
(overall  =  0.8,  ethylene/air  flames)  and  soot  sample  at  2.5  mm  (overall  =  0.8,  ethylene/air 
flames)  respectively.  The  coal  sample  was  15  %  lower  that  the  BET  value.  Two  soot  samples 
are  also  shown.  The  advantage  of  this  method  for  our  study,  as  compared  to  BET,  is  that  it 
allows  for  an  order  of  magnitude  less  sample.  Since  the  system  does  not  generate  a  lot  of  soot, 
this  was  an  important  consideration  for  these  experiments. 
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Figure  19.  Dubinnin-  Polanyi 
plots  showing  the 
experimental  data  plotted 
according  to  Equation  (18), 
and  the  linear  regressions  to 
obtain  the  monolayer 
adsorption  capacity  (xs)  for 
the  three  sample  run. 


Soot  Morphology  and  Nanostructure 

Samples  for  transmission  electron  microscopy  (TEM)  were  taken  using  a  thermophoretic  probe 
commonly  referred  to  as  a  “frog  tongue”  (see  Figure  20).  A  TEM  grid  holder  was  attached  to  a 
piston  and  compressed  air  at  60  psig  was  used  to  quickly  insert  the  TEM  grid  (200  mesh)  into  the 
flame  (Dobbins,  et  al.  1987).  Multiple  insertions  were  necessary  to  get  a  representative  soot 
sample  on  the  grid.  The  grid  was  oriented  with  the  face  perpendicular  to  the  gas  flow.  This  way 
allowed  for  sampling  at  the  exact  position  and  the  disturbance  of  the  flame  was  minimal.  Soot 
deposits  on  the  grid  because  of  the  thermophoretic  gradient  between  the  cold  grid  and  the  hot 
flame.  This  technique  allows  “freezing”  some  heterogeneous  reactions,  avoiding  changes  on  the 
soot  morphology  after  the  particles  have  impacted  upon  the  cold  surface.  HR-TEM  micrographs 
were  produced  using  two  FEI  TEMs,  Models  Tecnai  F30  and  F20  EFTEM  operated  under  200 
keV  accelerating  potential. 
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Figure  20.  Thermophoretic  sampling  system  for  transmission  electron  microscopy  analysis. 


HR-TEM  Image  Processing 

From  the  structural  data  obtainable  from  HR-TEM  images,  the  most  commonly  measured 
parameters  by  recent  studies  are  the  interlayer  distance,  d,  the  fringe  length,  L,  and  the  fractional 
coverage  of  detected  objects,  F  (Palotas  et  ah,  1996;  Sharma  et  ah,  1999;  Yang  et  ah,  2006). 
Except  for  F,  these  are  also  the  most  easily  verifiable  by  other  techniques,  such  as  X-Ray 
difffactometry  (Aso  et  ah,  2004).  For  the  direct  determination  of  d,  to  the  best  of  our  knowledge, 
all  authors  used  the  approximation  proposed  by  Palotas  et  ah  (1996),  where  the  interlayer 
distance  of  two  adjacent  carbon  fringes  is  the  vector- vector  distance  of  the  two  parallelized 
orientation  vectors.  These  vectors  are  obtained  by  averaging  their  orientations  and  keeping  their 
original  centers  of  mass.  Most  recent  publications  place  more  emphasis  on  the  classification  of 
adjacent  and  parallel  groups  of  fringes,  called  “stacks”.  Detailed  algorithms  for  deciding  an 
object’s  membership  to  stacks  are  presented  in  Sharma  et  ah  (1999)  and  Yang  et  ah  (2006). 
When  considering  data  describing  the  stacking  of  layers,  additional  parameters  arise  such  as  the 
diameters  of  stacks,  the  number  of  layers  in  stacks  and  other  derived  indices,  e.g.  the  crystallinity 
index  proposed  by  Yang  et  ah  (2006).  Using  whole  distribution  functions  of  d  and  the  fringe 
orientation  0  can  be  expressive,  because  the  amorphity  of  the  texture  -  the  deviation  from  perfect 
graphite  structure  -  can  be  described  by  the  statistical  indices  of  these  distributions.  Shim  and 
coworkers  (2000)  introduced  several  interesting  order  parameters,  such  as  the  2D  nematic  order 
parameter  S2,n  and  the  2D  polar  order  parameter  S2,p,  defined  as 


S  ZM  =  2cos2(0,)  -  1 

(21) 

$2,P  =  1  -  2cos2(fl,-) 

(22) 

Where  6  is  a  vector  containing  the  angles  between  the  orientation  vectors  of  each  fringe  and  a 
directional  vector;  the  overline  symbol  means  the  arithmetic  mean  value  of  the  elements  of  its 
argument.  The  directional  vector  is  the  reference  vector  for  S2;n  and  the  vector  pointing  from  the 
concentric  center  to  the  center  of  mass  of  fringe  i.  For  a  perfectly  ordered  graphite  structure, 
where  all  fringes  are  stacked  together  and  their  orientations  are  the  same,  the  nematic  order 
parameter  is  1 ;  for  a  perfectly  concentric,  onion-like  carbon  structure,  the  polar  order  parameter 
is  1. 

The  approach  used  in  this  study  is  fundamentally  different  from  the  previously  described  ones 
(Toth  et  al.  2010).  Instead  of  trying  to  separate  fringes  based  on  hypothesized  criteria  and  extract 
data  from  bitmaps,  we  tried  to  simplify  the  images  to  “frames”  that  can  be  topologically 
characterized.  These  topological  parameters  are  correlated  to  already  defined  properties.  For 
noise  filtration,  we  used  bandpass  gaussian  filters  in  the  frequency  domain.  Gaussian  filters  are 
better  than  ideal  filters  if  the  “ringing”  caused  by  the  sharp  changes  in  the  transfer  function  may 
lead  to  its  unwanted  detection,  as  in  our  case.  Despite  the  advantages  of  using  a  low-pass  filter  as 
proposed  by  Sharma  et  al.  (1999),  bandpass  filters  were  used,  because  of  the  additional  intensity 
homogenizing  properties  of  bandpass  Fourier  filtering  (inhomogeneous  luminance,  e.g.  large 
patches  of  darker  or  brighter  pixels  are  considered  very  low  frequency  noise).  The  calculated 
frequency  domains  were  centrally  shifted,  meaning  that  the  smallest  frequency  was  placed  in  the 
center  of  the  map.  Therefore,  the  coordinate  system  was  chosen  practically,  so  that  the  frequency 
of  a  point(xO;  yO)  of  the  frequency  map  was  calculated  using  the  following  equation: 
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where  x  and  y  are  the  vertical  and  horizontal  distances  of  a  point  from  the  center  of  the 
frequency  map.  Therefore,  a  two  dimensional  gaussian  transfer  function  with  a  mean  value  at  the 
central  point  is  defined  by  the  following  equation: 


0(x,y)  =  exp 
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In  the  case  of  a  Gaussian  distribution  function,  cr  would  be  the  variance.  Here  we  use  it  as  a 
parameter  controlling  the  suppression  of  unwanted  frequencies.  A  transfer  function  for  a 
bandpass  Gaussian  filter  can  be  generated  by  superimposing  two  Gaussian  filters,  with  the  high- 
pass  component  being  complemented.  Thus,  the  complete  function  can  be  defined  as  follows: 
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where  the  indices  HP  and  LP  mean  high-pass  and  low-pass,  respectively.  The  idea  behind  the 
determination  of  cr  is  that  the  same  extent  of  filtering  can  be  achieved  by  equal  “volume”  ideal 
and  Gaussian  transfer  functions.  The  ideal  bandpass  transfer  function  can  be  generally  pictured 
as  the  difference  of  two  cylinders;  their  radii  Dfip  and  Djmp  are  the  actual  high-  and  low-pass 
cut-off  frequencies  and  their  height  is  1.  Hence,  to  find  the  functions  a lp  =f  (DfLP)  and  Ohp  =f 
(D/hp)  for  I(x;  y),  an  mxn  size  image,  one  can  write 
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Considering  that  most  HR-TEM  images  are  digitized  in  mxm  sizes,  where  m  is  a  power  of  2  to 
ease  frequency  filtering,  so  m  =  n,  the  integration  of  the  left  sides  of  (26)  and  (27)  yields 
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It  is  apparent  that  the  two  equations  are  symmetrical,  therefore  both  <jlp  and  <yHp  can  be 
determined  with  the  implicit  transcendent  equation 
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where  XP  means  either  low-pass  (LP)  or  high-pass  (HP).  In  practical  situations,  when  Dfk  is 
significantly  lower  than  m,  Equation  (30)  can  be  simplified  to 
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The  relative  error  of  this  assumption  is  shown  in  Figure  21.  Note  that  the  maximum  frequency  of 
a  repeating  pattern  that  can  occur  on  an  mxm  size  image  is  V2m/2  and  the  values  of  the 
frequencies  that  need  to  be  filtered  increases  as  the  ratio  M  decreases.  In  our  case,  ctlp  and  ohp 
were  determined  at  cut-off  frequency  values  so  that  the  low-pass  component  filtered  out  patterns 
repeating  in  more  than  0.4  nanometers  and  the  high-pass  component  filtered  out  large  patches 
with  relative  frequencies  lower  than  10.  After  frequency  filtering,  the  inverse  Fourier 
transformed  image  was  re-scaled  and  saturated  to  eliminate  errors  caused  by  inhomogeneous 
illumination.  This  step  along  with  the  use  of  bandpass  Gaussian  transfer  functions  instead  of 
ideal  ones  was  necessary  because  of  the  characteristics  of  the  binarization  method.  For  the 
detection  process,  we  used  a  local,  adaptive  thresholding  technique.  Adaptive  local  threshold 
methods  have  been  widely  used  in  several  image  processing  areas,  e.g.  in  document  binarization 
and  handwriting  recognition  (Sauvola  et  al.,  2000;  Gatos  et  al.,  2006).  Local  means  that  the  value 
of  T  is  not  constant  for  the  entire  image;  instead,  the  image  is  divided  into  smaller  subregions, 
which  have  specific  values  of  T.  For  an  input  8  bit  grayscale  image  GR,  the  matrix  T  was 
calculated  in  the  following  way: 


£  e  \x  -  x  +fi]  D  (1,ml 
7j  e  [y  -  v, y  +  vl  n  [Ln] 

Gv  =  mcd  (g/^)  (32) 
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Figure  21.  The  relative  error  of  Equation  (31) 


Meaning  that  for  the  point  (x;  y),  tx,y  is  the  median  of  the  intensities  in  the  jUxv  neighborhood  on 
GR.  As  our  images  were  digitized  with  an  aspect  ratio  of  1 ,  square-shaped  median  filters  were 
used  (/d=  v).  From  our  point  of  view,  this  binarization  filter  that  produces  the  output  binary 
image  BW,  defined  in  the  following  way: 


Fu’vv 


0  ^  gJ.v.y  ^  fv.yi 

1  =>  grx.y  ^  G.y 


(33) 


There  are  several  advantages  to  this  method.  First,  median  filtering  has  noise-removal  effects, 
making  the  edges  of  the  detected  objects  smoother,  i.e.  it  tends  to  remove  spur  pixels,  which  is 
beneficial  in  the  next  step.  Second,  this  method  considers  regional  changes  in  intensity  rather 
than  absolute  values,  thus  less  focused  or  less  intense  fringes  can  be  detected  as  well.  Because  of 
this  feature,  Gaussian  frequency  filtering  and  re-saturation  are  crucial  steps  to  take  prior  to 
detection.  In  the  next  step,  the  detected  pattern  is  skeletonized  using  the  algorithm  presented  by 
Zhang  et  al,  (1984). This  iterative  algorithm  removes  edge  pixels  and  leaves  only  the  backbone  of 
the  structure,  but  does  not  tend  to  leave  spurs  and  diagonal  arms  (here  we  denote  this  operator  as 
“net”).  The  remaining  lines  make  a  binary  network  of  fringes  N  =  net  (BW),  which  is  the  input 
for  topological  analysis. 

Two  new  structural  parameters  are  proposed  that  can  be  measured  on  these  connected  networks. 
The  distance  deviation  parameter,  Q  [nm],  is  the  standard  deviation  of  the  Euclidean  distance  set 
E,  while  the  junction  parameter,  v  [1/nm],  is  the  average  number  of  branchpoints  per  network 
length  unit.  The  calculation  of  Q  can  be  achieved  in  a  number  of  steps.  First,  the  Euclidean 
distance  transform  of  the  binary  network  A  =  dist  (N)  is  computed.  The  Euclidean  distance 
transformed  value  of  a  point  bxy  in  a  binary  set  B  is  its  Euclidean  distance  to  the  nearest  point 
with  the  value  ~^bx;y.  Because  the  operator  ’net  ’  produces  binary  discrete  networks  without 
blocks  of  adjacent  1  values,  for  coordinates  (xo,yo),  where  nxo,yo=l,  Sxo:yo  is  always  a  regional 
minimum  of  A.  The  points  in  the  discrete  regional  maxima  in  A  correspond  to  half  of  the  shortest 
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distances  between  the  nearest  two  objects;  the  connected  lines  formed  by  the  adjacent  points  of 
the  regional  maxima  of  A  are  always  positioned  halfway  between  objects.  In  the  second  step, 
values  belonging  to  the  subset  of  regional  maxima  are  extracted  from  A.  Generally,  there  are  two 
obvious  ways  of  doing  this:  one  could  either  use  numerical  gradient  operators  to  find  regional 
maxima  or  apply  the  same  thinning  operator  ’net’  on  an  inverted  network  image.  The  use  of 
numerical  gradient  operators  was  found  to  provide  less  reliable  results,  as  it  tends  to  detect  pixels 
that  do  not  belong  to  the  above  mentioned  half  distances,  due  to  the  edge-enhancing  properties  of 

the  Laplacian  operator.  Instead,  the  coordinates  of  the  half-distances  were  found  by  the  thinning 

— ^ 

of  -,N.  The  Euclidean  distance  vector  E  is  defined  by  the  following: 

J  =  (fit.  A) :  net  (-.N)  =  1 1  nzn 


let  J *  be  a  vector  produced  by  reordering  the  elements  of  J  so  that 
./ ■  =  =  Afa  J  t  [1,p| 


(35) 


— ^ 

The  distance  deviation  parameter  is  defined  as  the  standard  deviation  of  the  elements  of  E  : 


iXm  ta-  -  E)2 
- 


(36) 


where  E  means  the  mean  value  of  the  elements  of  E .  As  can  be  seen,  Q  characterizes  the 
texture  in  a  complex  way,  as  it  is  a  parameter  describing  the  uniformity  of  the  lattice  spacing 
distribution,  the  parallelism  of  the  fringes  and  the  order  of  their  orientation  as  well.  A 
presentation  of  the  main  steps  of  the  algorithm  can  be  seen  in  Figure  22. 

To  present  the  correlation  of  Q  with  Sj,p  and  we  have  conducted  several  stochastic 
simulations  with  artificially  created  fringe  images.  For  the  analysis  of  Q  =  f  (SS.w),  the  images 
were  constructed  to  contain  short  fringes  with  randomly  positioned  centers  of  mass,  but  not 
allowing  their  contact.  Their  orientations  were  also  randomized,  by  using  a  random  generator 
yielding  uniformly  distributed  numbers  between  [-q;  q].  For  the  simulation  of  concentrically 
ordered  structures,  a  large  number  of  fringes  were  rendered  around  a  concentric  center,  equally 
spaced  from  each  other,  with  orientation  vectors  perpendicular  rotated  with  these  random  values. 
Examples  of  these  artificial  images  are  presented  in  Figure  23. 
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Figure  22. The  main  steps  of  the  new  algorithm  for  the  determination  of  distance  deviation 
parameter.  Left:  frequency  fdtered  image.  Center:  binarized  image  (BW)  with  the  thinned 
network  indicated  by  blue  lines  (N).  Right:  the  distance  transform  of  N,  A  and  the  extracted 
distances  indicated  by  blue  lines  (E).  The  padding  is  necessary  to  eliminate  errors  caused  by  the 

median  fdter  at  image  edges. 


Figure  23.  Examples  of  artificially  created  binary  fringe 
images.  Top  left:  oriented,  axially  symmetric  set.  S2;n~ 
0:99.  Top  right:  a  shued  version  of  the  oriented  axial  set. 
S2;N  ~  0:6.  Bottom  left:  oriented,  concentrically  symmetric 
set.  S2;p~  0:99.  Bottom  right:  a  shued  version  of  the 
oriented  polar  set.  S2;p~  0:6. 


If  we  simplify  Equation  (22)  to  only  consider  concentric  symmetry  and  the  number  of  fringes  on 
each  image  is  large,  both  and  Si,p  can  be  estimated  by 

SiX  *  E  {2  cos’  (0)  -  1) 

■  ;  (37) 

where  0  is  a  random  variable  of  uniform  distribution  with  values  of  [- q ;  q]  and  E  is  the  expected 
value  of  its  argument.  Solving  Equation  (37)  yields 


= 


sin(2r/) 

2? 


(38) 
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Figure  24  presents  the  results  of  the  stochastic  simulations.  For  each  parameter,  a  total  of  10  tests 
were  evaluated,  with  images  containing  400  fringes.  It  is  apparent,  that  the  decrease  in  the  order 
parameters  resulted  in  an  increase  in  Q,  but  the  correlation  is  not  linear.  The  measured  values  of 
Q  increased  monotonically  in  both  cases,  exhibiting  steeply  increasing  periods  at  both  ends  of 
the  curves.  With  our  sets,  changing  the  polar  order  parameter  of  the  concentrically  symmetric 
images  induced  less  predictable  changes  in  the  values  of  the  distance  deviation  parameter  and  the 
experimental  data  was  more  scattered.  Seemingly,  Q  responded  more  sensitively  to  the  changes 
of  the  axial  order  parameter.  The  monotonity  that  can  be  observed  in  Figure  24,  along  with  the 
fact  that  Q  responds  to  many  conditions  -  including  these  two  order  parameters  -  ensures  that 
even  slight  changes  in  structure  can  be  representatively  monitored.  The  downside  is  that  it  is 
impossible  to  gain  information  about  specific  structural  parameters  solely  by  observing  changes 
in  the  value  of  f>  it  should  rather  be  used  as  an  intuitive  empirical  parameter  of  orderliness. 

The  calculation  of  the  junction  parameter  is  obtained  by  counting  the  branchpoints  on  the  thinned 
networks  and  dividing  this  value  by  the  measured  length  of  the  network.  The  detected 
branchpoints  are  filtered  so  that  forks  with  an  arm  of  length  not  exceeding  a  certain  threshold  are 
not  considered  branches.  This  way  spur  pixels  and  arms  that  are  the  results  of  the  binarization  of 
jagged  contours  are  left  neglected.  The  filtering  algorithm  works  by  checking  each  branchpoints’ 
distance  to  the  nearest  endpoint  of  an  arm,  and  pixels  connecting  too  close  branch-end  pairs  are 
deleted  from  the  network.  Branchpoints  too  close  to  each  other  are  merged  by  binary  dilations. 

For  the  verification  of  the  junction  finding  algorithm,  we  used  HR-TEM  images  of  Si  110 
lattices.  At  certain  magnifications,  Si  110  single  crystals  look  like  perfect  grids,  where  intense 
spots  are  believed  to  be  Si  atoms.  For  such  arrangements  and  for  infinitely  large  fields  of  view, 
the  ideal  value  of  vcan  be  calculated  as  follows: 


v  =  lim 


m  ■  n 

n(m  -  1)L  +  m(n  -  1  )L 


(39) 


where  m  and  n  are  the  number  of  atoms  that  can  be  counted  on  the  image  horizontally  and 
vertically  and  L  ~  0:2715  nm  (O’mara,  1990)  (see  Figure  25). 
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Figure  24.  The  results  of  the  stochastic  simulations.  Top:  the  correlation  of  Q  with  S2N-  Bottom: 

the  correlation  of  with  S2;p. 
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Figure  25.  The  calibration  of  the  junction- finding  algorithm.  Left:  a  detail  of  a  HR-TEM 
micrograph  taken  of  an  oriented  Si  single  crystal.  Right:  The  detected  network  (red  lines)  and 

branchpoints  (green  dots). 


The  horizontal  axis  illustrates  the  size  of  the  detail  that  has  been  cut  from  a  large  lattice  image. 
Since  the  horizontal  and  vertical  dimensions  of  these  details  were  not  the  same,  an  average  size 
of 


m*  =  \'m  '  n 


(40) 


is  shown. 

Figure  26  shows  the  results  of  the  verification.  The  continuous  line  represents  the  ideal, 
calculated  value  of  v,  while  the  dots  mark  the  measured  results;  the  trend  can  be  easily  seen.  The 
errors  of  the  algorithm  are  most  commonly  caused  by  duplicate  branchpoints,  crystal 
imperfections,  inaccurate  network  extraction  and  image  noise,  however  the  magnitude  of  errors 
is  in  an  acceptable  range. 

Similar  to  distance  deviation,  v  is  an  empirical  parameter  of  the  orderliness  of  carbon  texture  as 
well:  higher  values  indicate  a  more  entangled  structure,  as  the  number  of  branchpoints  increase 
with  decreasing  fringe  lengths  and  overlapping  layers.  Its  value  is  0  for  an  ideal  graphitic  crystal 
with  parallel  fringes  and  no  overlapping. 
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Figure  26.  The  verification  of  the  junction-finding  algorithm  with  Si  1 10  lattice  images 


Methodology  to  predict  OH*  &  gas-phase  species  (H2,  O2,  CO,  CO2) 

Figure  27  represents  a  general  overview  of  the  methodology  to  predict  the  evolution  of  soot 
oxidizer  species  in  the  top  burner.  Experimental  measurements  of  major  gas-phase  species  (H2, 
CO,  CO2,  O2,  N2)  were  taken  within  the  top  burner.  The  results  of  these  measurements  at  the  top 
burner  surface  as  well  as  the  mass  flow  rate  and  the  temperature  profile  were  used  as  inputs  into 
CHEMKIN  4.1  to  predict  the  evolution  of  major  soot  oxidizers  (O2,  OH*)  in  the  top  burner.  The 
CHEMKIN  4.1  modeling  couples  a  lD-premixed  code  to  a  detailed  kinetic  model  defined  by  the 
user.  The  kinetic  mechanism  used  in  this  project,  developed  by  D’Anna  et  al.  (2006),  couples 
gas-phase  chemistry  with  a  sectional  approach  to  follow  the  transition  from  gas-phase  to  nascent 
particles  and  soot  oxidation.  The  details  of  the  model  have  been  reported  previously  (D’Anna  et 
al.,  2006,  2009).  Although  the  model  does  not  include  reactions  pathways  that  account  for  the 
fragmentation  of  particles,  its  gas-phase  chemistry  has  been  used  in  previous  investigations  to 
predict  the  concentration  profiles  of  radicals  such  as  H*,  O*,  OH*  etc.  in  premixed  flames. 

Methodology  to  Obtain  Soot  oxidation  Rates  from  Experimental  Data 

Soot  oxidation  rates  ( W)  were  calculated  from  the  experimental  PSDs  as  a  function  of  the  height 
above  the  second  burner  from: 


W 


8 


V  cm  sec  ) 


1  dm 
A  dt 


(41) 


where  A  was  the  surface  area  of  soot  particles  per  unit  volume  of  gas,  m  was  the  total  soot  mass 
per  unit  volume  of  gas,  and  t  was  time.  Area  and  mass  were  calculated  directly  from  the  SMPS 
diameter  and  number  concentration  assuming  spherical  particles.  Allowance  for  the  change  of 
soot  density  with  mobility  diameter  followed  the  procedures  of  Maricq  et  al.  (2004).  The  time 
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interval  was  obtained  from  measurements  of  temperature,  cross  sectional  area  and  volumetric 
flow  rate. 

These  data  were  used  to  predict  OH*  and  O2  concentrations  which,  in  turn,  were  used  in 
Equations  (1)  and  (2). 


Figure  27.  Schematic  representation  of  the  methodology  used  to  predict  the  evolution  of  soot 

oxidizer  species. 


RESULTS  AND  DISCUSSION 


Ethylene/air  Flames 

Temperature  profiles  were  measured  from  the  second  burner  surface  up  to  5  mm.  Figure  28 
illustrates  the  results  for  Ooverall  =  0.8  and  1.14.  The  difference  in  peak  temperature  is 
consistent  with  the  evolution  of  the  peak  temperature  in  premixed  flames  as  a  function  of  the 
equivalence  ratio.  These  temperature  profiles  were  used  as  inputs  into  a  detailed  kinetic  model  to 
predict  OH*  radical  concentration  profiles. 

Soot  size  distributions  were  obtained  as  a  function  of  the  HAB  for  the  lean  and  rich  flames 
studied.  Figure  29  presents  the  results  of  PSDs  for  Ooverall  =  0.8  (left)  and  Ooverall  =1.14 
(right)  conditions  in  the  locations  where  major  changes  took  place.  Data  for  the  leanest  flame 
(Ooverall  =  0.8)  at  HAB  =  0  mm  (Figure  29,  left),  showed  that  the  initial  distribution  was  mostly 
in  the  fine  mode  (particles  diameter  >10  nm)  with  little  contribution  of  particles  less  than  10  nm. 
This  PSD  did  not  change  significantly  up  to  2.5  mm  where  particles  in  the  ultrafine  mode 
(particles  diameter  <10  nm)  started  to  appear.  In  the  range  of  HAB  from  2.5  to  3  mm,  particle 
mean  diameter  decreased  from  36.7  nm  to  3.7  nm,  and  at  higher  HAB,  the  particles  started  to 
burnout  significantly,  evidenced  by  a  drop  in  number  concentration.  For  Ooverall  =1.14  (Figure 
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29,  right),  particles  in  the  ultrafine  mode  almost  disappeared  when  they  reached  HAB  =  2  mm, 
and  they  are  difficult  to  see  at  higher  elevations.  For  the  fine  mode,  the  rate  of  oxidation 
increased  significantly  after  HAB  =  2  mm,  and  soot  burnout  dominated  the  oxidation  process. 
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Figure  28.  Top  burner  temperature  profiles  as  a  function  of  HAB  for  Ooverall  =  0.8  and  1.14. 

Thermocouple  uncertainty  ±  5  OK 

Particles  in  size  intervals  of  less  than  and  greater  than  10  nm  (the  detection  limit  of  the  SMPS 
was  3  nm)  were  integrated  to  provide  measures  of  the  number  and  mass  concentration  of  the 
ultrafine  and  fine  modes  (Figure  30a  and  30b,  for  the  lean  flame  and  Figure  30c  and  30d  for  the 
rich  flame).  For  the  lean  flame,  the  particle  number  (Figure  30a)  and  mass  concentration  (Figure 
30b)  in  the  ultrafine  mode  were  observed  to  increase  up  to  HAB  =  3  mm.  This  behavior  has  been 
attributed  to  oxygen  diffusion  into  the  pores  of  the  soot  particles  causing  internal  burning,  which, 
in  turn  caused  the  soot  aggregates  to  break  apart  (Neoh  et  al.,  1981,  1985;  Fenimore  et  al.,  1967). 
It  is  also  interesting  to  notice  that  this  fragmentation  took  place  for  soot  burnout  less  than  10%, 
which  is  considerably  lower  than  the  value  observed  by  Neoh  and  coworkers  (80%  burnout).  At 
HAB  >  3  mm,  particles  in  the  ultrafine  mode  were  almost  completely  combusted,  as  evidenced 
by  the  decrease  in  number  and  mass  concentration.  For  the  fine  mode  (Figure  30a  and  30b),  mass 
and  number  remained  relatively  constant  up  to  2.5  mm,  and  then  started  to  decrease  in  the  region 
between  2.5  to  3  mm;  this  can  be  attributed  to  the  fragmentation  and  oxidation  of  the  bigger 
particles.  In  the  upper  region  of  the  flame,  soot  was  almost  completely  burned  out.  For  the 
richest  flame  (Figure  30c  and  30d),  both  ultrafine  and  fine  modes  showed  a  decrease  in  number 
and  mass  concentration  up  to  3  mm.  Higher  in  the  burner,  they  leveled  off  with  only  minor 
changes  in  number  or  mass  concentration.  No  evidence  of  fragmentation  was  observed  for  this 
flame  in  the  range  of  HAB  from  0  to  5  mm. 
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Figure  29.  Experimental  PSDs  measured  as  a  function  of  the  height  above  the  top  burner  for 
Ooverall  =  0.8  (left)  and  Ooverall  =  1.14  (right). 


39 


Figure  30.  Integrated  number  and  mass  concentrations  for  fuel  lean  and  fuel  rich  flames, 
a)  Number  concentration  as  a  function  of  HAB  for  the  ultrafine  and  fine  modes  (d>0Veraii 
=  0.8).  b)  Mass  concentration  as  a  function  of  HAB  for  the  ultrafine  and  fine  modes 
(d1  overall  =  0.8).  c)  Number  concentration  as  a  function  of  HAB  for  the  ultrafine  and 
fine  modes  (d>0Veraii  =  1.14).  d)  Mass  concentration  as  a  function  of  HAB  for  the 
ultrafine  and  fine  modes  (0OVeraii  =  1.14). 

The  concentration  profiles  of  major  gas  species,  H2,  O2,  CO  and  CO2,  were  measured 
experimentally  within  the  secondary  flame  using  the  isokinetic  sampling  system  coupled  to  a 
gas-chromatograph  described  in  the  Materials  and  Methods  section.  Experimental  temperature 
profiles  (Figure  28),  and  experimental  conditions  measured  at  the  surface  of  the  top  burner  (mass 
flow  rate,  H2,  CO,  CO2,  O2  concentration,  and  PSDs)  were  used  as  inputs  into  the  detailed  kinetic 
model  (D’Anna  et  al.)  to  predict  the  OH*  concentrations  as  a  function  of  the  HAB.  Figure  31a 
and  31b  compare  the  experimental  measurements  (open  symbols)  of  H2,  O2,  CO  and  CO2  to 
model  predictions  (solid  and  dashed  line)  for  ^overall  of  0.8  and  1.14  respectively.  Predictions 
of  OH*  are  also  presented  and  they  are  plotted  with  the  right-side  scale.  Good  agreement  was 
obtained  between  the  model  results  and  the  experimental  data  for  these  flames.  Sampling  effects 
due  to  the  water-cooled  probe  were  taken  into  account  by  shifting  the  model  downstream  as 
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supported  in  previous  investigations  (Echavarria  et  al.,  2009;  Zhao  et  al.,  2003).  As  seen  in  these 
figures,  at  the  burner  surface  (HAB  =  0  mm),  the  main  fuel  components  for  the  secondary, 
premixed  flame  were  H2  and  CO.  In  the  leanest  flame,  Ef  and  CO  were  consumed  in  the  upper 
region  of  the  flame,  oxygen  dropped  to  about  6  %,  and  CO2  concentrations  increased  and 
leveled.  According  to  predictions,  OH*  started  to  increase  in  the  region  where  O2  decreased  and 
reached  a  maximum  when  the  concentration  of  O2  leveled  off  for  HAB  >  4  mm.  The 
concentration  of  O2,  CO  and  CO2  started  to  show  a  significant  change  after  an  induction  period 
up  to  2.5  -  3  mm  which  is  consistent  with  the  trend  of  the  PSDs  obtained  in  the  second  burner, 
where  major  changes  in  number  and  mass  concentration  occurred.  On  the  other  hand,  H2,  CO, 
O2,  CO2  and  OH*  were  consumed  and  formed  closer  to  the  burner  surface  for  the  rich  condition; 
the  higher  OH*  concentrations  favor  faster  soot  oxidation  resulting  in  little  or  no  fragmentation 
of  the  particles  in  this  flame.  The  effect  of  other  radicals  and  species,  such  as  O*,  H2O,  CO2,  and 
CO,  etc.  on  soot  oxidation  were  not  considered  since  previous  studies  (Neoh  et  al.,  1981,  1985; 
Vierbaum  et  al.,  2003)  have  shown  that  these  species  contributed  negligibly  to  the  oxidation 
under  experimental  conditions  similar  to  those  in  the  present  experiments. 
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Figure  3 1 .  Experimental  measurements  versus  model  predictions  of  H2  (□  data,  —  model),  O2 

(A  data, - model),  CO  (o  data, - model),  CO2  (x  data,  —  model)  and  OH*(. - model) 

for  (a)  Ooverall  =  0.8  and  (b)  Ooverall  =  1.14.  Error  bars  represent  the  experimental  uncertainty. 


The  results  for  the  experimental  W  calculated  using  Equation  (41)  are  presented  in  Figure  32  for 
Ooverall  =  0.8  and  1.14.  For  Ooverall  =  0.8,  rates  were  found  to  be  relatively  low  and  constant 
close  to  the  burner  surface,  increasing  in  the  upper  (HAB  >  2.5  mm)  region  of  the  flame.  For  the 
rich  flame  (Ooverall  =  1.14),  rates  increased  faster  close  to  the  burner  surface  and  reached  a 
maximum  of  7.8xl0~4  g/cm2-sec  before  starting  to  decrease.  This  behavior  was  consistent  with 
the  trend  of  the  PSDs  (Figure  29),  and  O2  and  OH*(see  Figure  31)  concentrations  profdes  where 
small  changes  took  place  up  to  1 .5  -2  mm. 
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Figure  32.  Experimental  soot  oxidation  rates  (W)  calculated  from  the  evolution  of  PSD  for 

Ooverall  =  0.8  and  1.14. 


Soot  generated  in  the  first  burner  in  a  ethylene/air  flame  was  oxidized  not  only  for  the  conditions 
presented  above  (<E>overall=  0.8  and  1.14)  but  for  intermediate  conditions  that  allowed  further 
study  of  the  effect  of  temperature  and  equivalence  ratio  on  soot  oxidation.  Similarly  to  the 
previous  data,  temperature,  PSDs,  number  and  mass  concentration,  gas-phase  compounds  and 
soot  oxidation  rates  were  characterized  in  the  two-stage  burner  as  a  function  of  the  HAB  for 
Ooverall  =  0.87,  0.94,  and  1.07. 

Figure  33  illustrates  the  temperature  profile  for  these  experimental  conditions  measured  within 
the  second  burner  up  to  5  mm.  Temperature  profiles  showed  a  typical  behavior  of  the  evolution 
of  peak  temperature  with  equivalence  ratio.  Peak  temperature  for  the  lean  flames  increased  from 
the  leanest  conditions  0.8  (see  Figure  28)  to  0.94  (Figure  33).  On  the  other  hand,  temperature 
decreased  in  the  rich  side  from  Ooverall  =  1.07  to  1.14.  Once  again,  these  temperature  profiles 
were  used  as  inputs  into  the  detailed  kinetic  model  to  obtain  the  concentration  of  OH*  as  a 
function  of  the  HAB. 


Figures  34a,  34b  and  34c  present  the  results  for  the  evolution  of  the  soot  size  within  the  flame 
from  0  to  5  mm  for  Ooverall  =  0.87,  0.94  and  1.07  respectively.  These  PSDs  were  integrated  in 
size  intervals  of  less  than  and  greater  than  10  nm  to  characterize  the  number  and  mass 
concentrations  in  the  ultrafine  (3-10  nm)  and  fine  (>  10  nm-160)  modes  (Figure  35a,  b,  c  and  d 
for  Ooverall  =  0.87,  0.94  and  Figure  35e  and  f  for  Ooverall  =  1.07).  Fragmentation,  to  a  smaller 
extent,  was  only  observed  to  occur  for  the  0.87  flame.  However,  it  seemed  to  disappear  for  the 
0.94  and  1.07  flames  where  soot  burnout  dominated  the  oxidation  process.  These  results  were 
reflected  in  terms  of  number  and  mass  concentration  in  both  the  ultrafine  and  fine  modes  (Figure 
35).  Number  and  mass  concentration  increased  close  to  the  burner  surface  (HAB  <1.5  mm)  for 
Ooverall  =  0.87  in  the  ultrafine  mode,  while  higher  in  the  burner,  both  number  and  mass 
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concentration  showed  a  continuous  decrease,  and  started  to  level  after  6  mm.  In  the  fine  mode, 
close  to  the  burner  surface  (HAB  <1.5  mm),  a  slight  decrease  in  mass  was  observed  which  was 
consistent  with  the  increase  in  mass  in  the  ultrafine  mode.  Higher  in  the  burner,  the  soot  burnout 
process  dominated  in  both  the  ultrafine  and  fine  modes.  For  equivalence  ratios  of  0.94  and  1.07 
(Figure  35c,  d,  e  and  f),  both  the  ultrafine  and  fine  modes  showed  major  changes  in  mass  and 
number  concentration  after  1.5  mm  where  the  rate  of  oxidation  increased  significantly  (see 
Figure  35a),  probably  as  a  result  of  increased  surface  area  with  the  smaller  particles  or  other 
resulting  structural  changes. 


1600 


c\ 

U 

3 

u 

QJ 

a 

S 

QJ 

H 

QJ 

s 


1500 


1400 


1300 


1200 


1100 


1000 


L_J 


:|f;3s  s  □ 
*  $-•&  & 


O 


overall 


1.07 


;3§t-.  s  s  s 


*  0.87 
0  0.94 


1  2  3  4  5 

HAB,  mm 

Figure  33.  Top  burner  temperature  profiles  as  a  function  of  HAB  for  ^overall  =  0.87,  0.94  and 

1.07.  Thermocouple  uncertainty  ±  50K 


In  terms  of  the  gas-phase  compositions  (Figure  36),  the  evolution  of  the  species  such  as  H2,  O2, 
CO,  CO2  and  OH*  was  similar  to  the  results  presented  for  ^overall  =  0.8,  and  1.14.  Once  again, 
it  should  be  highlighted  the  earlier  formation  of  the  OH*  in  the  rich  flame  (1.07),  which 
promoted  the  burnout  process  instead  of  fragmentation.  Although  fragmentation  may  be  taking 
place  close  to  the  burner  surface,  the  higher  reactivity  of  the  OH*  favors  the  faster  consumption 
of  the  soot  particles. 

Soot  oxidation  rates  (W)  (see  Figure  37)  were  also  calculated  from  the  evolution  of  PSDs.  For 
^overall  =  0.87,  0.94,  and  1.07  rates  were  higher  closer  to  the  burner  surface  relative  to  the  Ws 
estimated  for  the  leanest  flame  (Ooverall  =  0.8)  (see  Figure  32).  The  little  or  no  fragmentation 
found  in  the  flames  0.87,  0.94  and  1.07  can  be  supported  for  the  observed  increase  in  soot 
oxidation  rates  and  the  relatively  high  concentration  of  OH*  higher  in  the  burner.  As  mentioned 
before,  the  higher  reactivity  of  the  OH*  promoted  the  earlier  consumption  of  particles  which  was 
reflected  in  terms  of  the  evolution  of  PSDs,  number  and  mass  concentration  within  the  flames. 


43 


1  10  100  1000 


Particle  Diameter,  nm 


Figure  34.  Experimental  PSDs  measured  as  a  function  of  the  height  above  the  top  burner  for  (a) 
Ooverall  =  0.87,  (b)  Ooverall  =  0.94  and  (c)  Ooverall  =  1.07. 
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Figure  35.  Integrated  number  and  mass  concentrations  for  fuel  lean  and  fuel  rich  flames. 
(a)and  (b)  number  and  mass  concentration  for  the  ultrafine  and  fine  modes  (Ooverall  =  0.87). 
(c)  and  (d)  number  and  mass  concentration  for  the  ultrafine  and  fine  modes  (Ooverall  =  0.94). 
(e)  and  (f)  number  and  mass  concentration  for  the  ultrafine  and  fine  modes  (Ooverall  =  1.07). 
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Figure  36.  Experimental  measurements  (open  symbols)  versus  model  predictions  of  FE  (□,  data, 

—  model),  O2  (A  data, - model),  CO  (o  data, - model),  CO2  (X  data,  —  model)  and 

OH*(. - model)  for  (a)  ^overall  =  0.87,  (b)  ^overall  =  0.94  and  (c)  ^overall  =  1.07.  Error 

bars  represent  the  experimental  uncertainty. 
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Figure  37.  Experimental  soot  oxidation  rates  (W)  calculated  from  the  evolution  of  PSD  as  a 
function  of  the  HAB  for  (Doverall  =  0.87,  0.94  and  1.07. 


Surrogate/ Air,  m-Xylene/Air  and  n-Dodecane/Air  Flames 

The  use  of  surrogate/air,  m-xylene/air  and  n-dodecane/air  flames  constituted  our  transition  from 
simple  fuels,  i.e.,  ethylene,  to  more  practical  liquid  fuels.  Preliminary  studies  with  a  mixture 
23/77  vol.  %  were  difficult  to  perform  due  to  the  high  sooting  tendency  of  the  23/77  mixture. 
Studies  on  soot  oxidation  were  also  carried  out  for  the  pure  components  of  the  surrogate  mixture 
m-xylene  and  n-dodecane. 

The  oxidation  burner  was  initially  characterized  in  terms  of  temperature  under  lean  (^overall  = 
0.8)  and  rich  (Ooverall  =  1.14)  conditions.  Flame  temperature  profiles  measured  in  the  top 
burner  for  conditions  where  soot  was  derived  from  m-xylene/air,  surrogate/air,  and  n- 
dodecane/air  flames  are  presented  in  Figures  38a, b  and  c,  respectively.  Temperature  profiles  for 
the  six  flames  showed  a  similar  behavior,  increasing  rapidly  close  to  the  burner  surface  until  a 
peak  temperature  was  reached.  Downstream,  they  leveled  off  and  remained  relatively  constant  up 
to  7  mm,  which  is  the  maximum  distance  above  the  second  burner  surface  where  experimental 
measurements  were  taken.  Temperature  profiles  and  experimental  conditions  at  the  top  burner 
surface  were  used  as  inputs  into  CHEMKIN  4.1  to  model  the  evolution  of  major  gas-phase 
species  (H2,  O2,  CO,  CO2,  OH*). 
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Figure  38.  Top  burner  temperature  profiles  as  a  function  of  HAB  under  lean  and  rich  conditions. 

Soot  in  the  top  burner  was  derived  from:  (a)  m-xylene/air  flame,  (b)  surrogate/air  and  (c)  n- 
dodecane/air  flames.  Error  bars  represent  thermocouple  uncertainty  (±  50K) 

The  oxidation  process  of  soot  particles  formed  in  the  bottom  burner  was  mainly  followed  by 
measuring  the  evolution  of  particle  size  distribution  (PSD)  and  number  concentration  as  a 
function  of  the  HAB.  Figure  39a  (0  to  3.5  mm)  and  Figure  39b  (4  to  7  mm)  illustrates  the  results 
for  the  lean  (Ooverall  =  0.8)  surrogate/air  flame  where  the  most  significant  changes  in  terms  of 
PSD  took  place.  At  the  top  burner  surface  (HAB  =  0  mm),  the  PSD  was  unimodal  and  most  of 
the  distribution  was  dominated  by  the  fine  mode  (particle  diameter  >10  nm).  In  this  flame,  the 
PSD  evolved  from  a  unimodal  to  bimodal  distribution  with  the  smaller  nanoparticles 
corresponding  to  the  ultrafine  mode  (particle  diameter  <  10  nm).  The  presence  of  these 
nanoparticles  in  the  ultrafine  mode  was  mainly  associated  with  the  fragmentation  of  the  bigger 
particles  at  low  burnout  (~  15  %).  As  mentioned  before,  previous  studies  (Neoh,  1981;  Merrill, 
2005;  Lighty,  2008;  Echavarria,  2011)  have  attributed  this  behavior  to  O2  penetrating  the  soot 
nanostructure,  which  in  turn,  caused  internal  burning  and  led  to  the  fracture  and  breakup  of  the 
particles.  In  the  upper  regions  of  the  flame  (HAB  >  3.5  mm,  Figure  39b),  number  concentration 
decreased  in  both  the  ultrafine  and  fine  mode,  and  at  6.5  mm  above  burner  surface  most  of  the 
particles  were  burned  out.  On  the  other  hand,  the  series  of  oxidation  experiments  performed  on 
soot  derived  from  the  surrogate  flame  under  rich  (Ooverall  =  1.14,  Figures  39c  and  39d) 
conditions  also  showed  the  presence  of  particle  fragmentation,  a  characteristic  that  was  not 
observed  in  previous  studies  when  soot  derived  from  methane/air  (Neoh,  1981)  or  ethylene/air 
(Echavarria  et  al.,  201 1)  flames  was  oxidized  under  slightly  rich  conditions. 
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Figure  39.  Evolution  of  PSD  in  the  top  burner  in  the  surrogate/air  flame  for  (a)  ^overall  =  0.8 
and  HAB  from  0  to  3.5  mm,  (b)  ^overall  =  0.8  and  HAB  from  4  to  7  mm,  (c)  Ooverall  =1.14 
and  HAB  from  0  to  4  mm  and  (d)  ^overall  =1.14  and  HAB  from  5  to  7  mm 


Results  in  terms  of  the  evolution  of  soot  size  distribution  in  the  top  burner  under  lean  and  slightly 
rich  conditions  for  soot  derived  from  rich  flames  of  the  pure  components  of  the  surrogate 
mixture  are  presented  in  Figures  40a  through  40g.  Fragmentation  occurred  in  all  the  cases  as 
reflected  for  the  presence  of  nanoparticles  in  the  ultrafine  mode  after  a  characteristic  initial 
induction  period  between  1  to  2  mm.  Higher  in  the  burner  most  of  the  particles  were  burnout 
under  fuel  lean  conditions  and  the  concentration  of  the  fine  sized  particles  remains  relatively 
high  for  rich  conditions. 

These  results  were  more  clearly  seen  when  the  PSDs  were  integrated  in  size  intervals  of  less  than 
and  greater  than  10  nm  to  provide  measures  of  the  change  in  number  concentration  in  the 
ultrafine  and  fine  modes  as  a  function  of  the  HAB  (Figures  41a  and  b,  lean  and  rich  ultrafine 
mode  and  41c  and  d  lean  and  rich  fine  mode).  Figures  41a  and  b  clearly  shows  the  drastic 
increase  of  nanoparticles  in  the  ultrafine  mode  close  to  the  burner  surface  for  lean  and  rich 
conditions,  while  higher  in  the  burner,  number  concentration  in  both  the  ultrafine  and  fine  modes 
(Figures  41c  and  d)  showed  a  continuous  decrease,  particularly  for  the  lean  flames.  For  the  rich 
flames  (Figure  41b),  it  is  clear  that  the  extent  of  fragmentation,  which  is  characterized  by  the 
increase  of  nanoparticles  in  the  ultrafine  mode,  was  higher  relative  to  lean  (Figure  41a) 
conditions,  and  it  lasted  over  almost  the  entire  range  of  HABs  evaluated.  The  decrease  of 
nanoparticles  in  the  fine  mode  under  the  rich  condition  was  also  observed  to  be  considerably 
lower  relative  to  lean  flame  (Figure  41c  and  d).  In  addition,  the  figures  show  that  the  surrogate 
does  not  appear  to  burnout  as  quickly  as  the  pure  components,  especially  in  the  lean  case. 
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Figure  40.  Evolution  of  PSD  in  the  top  burner  in  the  systems  m-xylene/H2/air  and  n- 
dodecane/H2/air  for  (a)  m-x/FE/air,  Ooverall  =  0.8  and  HAB  from  0  to  3  mm,  (b)  m- 
x/H2/air  ,  Ooverall  =  0.8  and  HAB  from  4  to  5.5  mm,  (c)  m-x/H2/air,  Ooverall  =1.14  and  HAB 
from  0  to  3  mm,  (d)  m-x/H2/air  ,  Ooverall  =1.14  and  HAB  from  3.5  to  6  mm,  (e)  n-d/H2/air, 
Ooverall  =  0.8  and  HAB  from  0  to  3  mm,  (f)  n-d/H2/air  , Ooverall  =  0.8  and  HAB  from  3.5  to  6 
mm,  (g)  n-d/H2/air,  Ooverall  =1.14  and  HAB  from  0  to  3  mm,  (h)  n-d/H2/air ,  Ooverall  =1.14 

and  HAB  from  4  to  7  mm. 
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Figure  41.  Number  concentration  measurements  in  the  ultrafine  (3-10nm)  and  fine  (>10nm) 
modes  calculated  from  the  PSD  for  m-xylene,  surrogate,  and  n-dodecane  flames  under  fuel  lean 
and  fuel  rich  condition,  (a)  ^overall  =  0.8,  ultrafine  mode,  (b)  Ooverall  =  1.14,  ultrafine  mode, 
(c)  Ooverall  =  0.8,  fine  mode  and  (d)  Ooverall  =  1.14,  fine  mode.  *  Ooverall  under  fuel  lean 
conditions  for  the  system  with  m-xylene/tb/air  is  equal  to  0.87. 

The  mole  fractions  of  gas-phase  species  were  measured  using  online  GC  analysis.  OH*  was 
calculated  by  CHEMKIN  4. 1  software  as  previously  discussed.  However,  in  this  case  a  validated 
JP-8  mechanism  developed  by  Zhang  et  al.  (2005)  was  used  to  predict  the  exit  stream  from  the 
first  burner.  The  results  showed  that  the  most  important  species  that  mixed  with  the  secondary  air 
were  H2,  CO,  CO2,  H2O,  and  soot.  The  second  burner  was  modeled  using  a  detailed  kinetic 
model  developed  by  D’Anna  et  al.  (2006,  2008)  that  includes  reaction  paths  for  the  oxidation  of 
H2,  CO  and  the  formation/oxidation  of  soot  in  flames.  Experimental  measurements  of  gas  flow 
rates,  concentration  of  gas-phase  species  (H2,  CO,  CO2,  H2O,  O2,  soot)  at  the  burner  surface  and 
the  temperature  profile  were  used  in  the  simulation  of  the  top  burner.  Figures  42a  through  42i 
present  the  measured  and  predicted  profiles  of  H2,  CO,  CO2,  O2  and  OH*  within  the  top  burner 
for  lean  (Ooverall  =  0.8,  0.87)  and  rich  (Ooverall  =  1.14)  conditions  in  the  surrogate/air  (Figure 
42a,  b  and  c),  m-xylene/air  (Figure  42d,  e  and  f)  and  n-dodecane/air  (Figure  42g,  h  and  i) 
systems.  Model  results  were  shifted  downstream  to  account  for  sampling  probe  effects  (Zhao  et 
al.,  2003,  Echavarria  et  al.,  2009).  Among  the  main  observations,  in  the  surrogate  mixture 
system  (Figures  42a,  lean  and  42b,  rich)  the  concentration  profiles  of  H2,  CO,  CO2  and 
particularly  O2  showed  a  major  change  after  an  induction  period  up  to  3  mm  and  4  mm  for  lean 
and  rich  flames  respectively.  This  behavior  was  consistent  with  the  regions  where  the  most 
significant  changes  in  terms  of  size  and  number  concentration  in  both  flames  were  observed. 
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The  high  temperature  environment  and  the  relatively  high  concentration  of  O2  in  that  region 
favored  the  fragmentation  of  soot  particles  for  the  lean  flame.  It  is  also  important  to  note  that  for 
the  rich  condition,  a  more  prolonged  induction  period  took  place,  which  in  turn,  favored  the 
higher  extent  of  fragmentation  under  this  condition.  Figure  42c  shows  the  calculated 
concentration  profiles  of  OH*  radicals  within  the  top  burner  for  the  same  system.  The  evolution 
of  these  species  was  similar  for  both  lean  and  rich  flames.  The  concentration  of  OH*  had  a  major 
impact  in  the  upper  regions  of  the  flames  where  it  contributed  mostly  to  the  burnout  of  particles, 
in  particular  for  the  lean  condition  where  the  concentration  of  OH*  was  relatively  higher.  The 
major  drop  of  particle  concentration  in  the  upper  regions  of  the  flames  (see  Figure  41)  was 
attributed  to  the  higher  intrinsic  reactivity  of  OH*  compared  to  O2  (Neoh,  1981).  For  the  other 
two  systems  where  soot  was  derived  from  m-xylene/air  (Figures  42d  through  42f)  and  n- 
dodecane/air  (Figures  42g  through  i)  flames,  the  evolution  of  major  gas-phase  species, 
particularly  the  changes  in  O2  and  OH*  showed  a  similar  behavior.  O2  and  OH*  changed  in  the 
region  where  the  most  significant  changes  in  terms  of  PSD  and  number  concentration  took  place. 
O2  was  predominant  on  soot  oxidation  close  to  the  burner  surface  and  OH*  in  the  upper  regions 
of  the  flames.  In  addition,  the  reduced  amount  of  OH*  predicted  in  the  upper  regions  of  the 
flame  could  account  for  the  reduced  rate  of  fine  particle  decrease  for  the  surrogate  as  compared 
to  m-xylene  and  dodecane. 

Soot  oxidation  rates  (W)  were  calculated  from  the  evolution  of  PSDs  using  Equation  41  (see 
Figure  43).  For  lean  and  rich  conditions,  the  rates  were  relatively  higher  than  the  soot  oxidation 
rates  calculated  for  ethylene/air  flames.  However  they  are  still  lower  closer  to  the  burner  surface 
and  showed  a  continue  increase  with  HAB.  The  higher  soot  oxidation  rates  obtained  upper  in  the 
burner  were  associated  to  the  higher  reactivity  of  the  OH*  Neoh  et  al.  1981)  that  promoted  the 
consumption  of  particles  higher  in  the  burner  for  lean  and  rich  conditions. 
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Figure  42.  Evolution  of  major  gas  phase  species  for  the  surrogate  mixture  and  the  surrogate 
components  within  the  second  burner.  Symbols  correspond  to  experimental  results  and  solid  and 
dashed  lines  to  model  predictions  (a)  H2,  O2,  CO  and  CO2  for  surrogate,  (b)  H2,  O2,  CO  and  CO2 
for  surrogate,  (c)  OH*  predictions,  surrogate,  Ooverall  =  0.8  and  1.14;  (d)  H2,  O2,  CO  and  CO2 
for  m-xylene,  (e)  H2,  O2,  CO  and  CO2  for  m-xylene,  (f)  OH*  predictions,  m-xylene,  Ooverall  = 
0.87  and  1.14;  (g)  H2,  O2,  CO  and  CO2  for  n-dodecane,  (h)  H2,  O2,  CO  and  CO2  for  n-dodecane, 
(i)  OH*  predictions,  n-dodecane,  Ooverall  =  0.8  and  1.14.  Error  bars  represent  experimental 

uncertainty 
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Figure  43.  Experimental  soot  oxidation  rates  (W)  calculated  from  the  evolution  of  PSD  as  a 
function  of  the  HAB  under  lean  (<|)=0.8)  and  rich  (4>=1.14)  conditions. 


Development  of  kinetics  for  soot  oxidation  and  comparisons 

Experimental  soot  oxidation  rates  estimated  from  SMPS  measurements  and  presented  in  Figures 
32,  37  (for  ethylene/air)  and  43  (for  the  surrogate  mixture  and  its  pure  components)  were 
compared  to  previously-mentioned  models  used  in  the  literature  to  account  for  soot  oxidation  via 
O2  (NSC,  Equation  2)  and  soot  oxidation  via  OH*  (Neoh,  Equation  1).  The  data  are  presented 
with  the  temperature  profiles  corrected  for  probe  effects  and  O2  and  OH*  concentrations 
corresponding  to  these  temperatures.  Each  is  shown  separately  but  the  entire  oxidation  is  the 
addition  of  the  two  effects.  The  results  under  lean  and  rich  overall  conditions  are  presented  in 
Figures  44a  through  44h  for  the  combustion  systems  where  soot  was  derived  from  ethylene/air 
(Figures  44a  and  44b),  surrogate/air  (Figures  44c  and  44d),  m-xylene/air  (44e  and  44f)  and  n- 
dodecane/air  (Figures  44g  and  44h)  premixed  flames. 

Figure  44  shows  that  under  the  experimental  conditions  used  in  this  work,  the  NSC  equation 
under-predicted  the  set  of  experimental  data  in  the  region  where  soot  oxidation  via  O2  became 
predominant  as  evidenced  by  the  gas-phase  analysis  (see  Figures  31  and  42)  for  the  lean  cases. 
In  the  fuel-rich  cases,  the  predictions  were  much  lower.  Neoh’s  kinetic  equation  to  account  for 
soot  oxidation  via  OH*  produced  a  relatively  good  fit  higher  in  the  burner,  which  confirmed  that 
the  most  significant  contribution  of  OH*  radical  to  soot  burnout  was  in  that  region.  It  should  be 
noted  that  the  error  in  temperature  is  likely  the  highest  uncertainty  and  HAB  could  be  +/-  1mm. 
Under  lean  conditions,  oxidation  via  O2  is  also  occurring  in  the  upper  regions  of  the  flame; 
however,  as  expected,  because  of  the  higher  intrinsic  reactivity  of  OH*  with  soot  (Neoh,  1981), 
oxidation  via  O2  had  a  minimum  impact  in  the  upper  regions  of  the  flame.  Soot  oxidation  via  O2 
contributed  mainly  to  soot  fragmentation  in  the  regions  close  to  the  burner  surface  and  higher  in 
the  top  burner  for  the  lean  oxidation  flames. 
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Figure  44.  Experimental  versus  models  prediction  using  NSC  and  Neoh’s  kinetic  oxidation 
equation  for  the  set  of  flames  ran  under  lean  and  slightly  rich  conditions,  (a)  and  (b)  ethylene/air; 
(c)  and  (d)  surrogate/air;  (e)  and  (f)  m-xylene/air  and  (g)  and  (h)  n-dodecane/air. 
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Comprehensive  Soot  Oxidation  Expression 

Neoh’s  (1981)  kinetic  expression  to  account  for  soot  oxidation  via  OH*  has  been  used 
extensively  in  the  literature  and  in  this  work  as  well.  The  results  presented  in  the  previous  section 
showed  that  this  expression  was  suitable  to  account  for  the  oxidation  of  particles  derived  from 
different  sources  under  fuel  lean  and  slightly  rich  conditions  when  oxidation  via  OH*  was 
dominant.  A  kinetic  expression  that  includes  the  effect  of  temperature,  O2  and  OH* 
concentration  to  calculate  the  oxidation  rate  of  soot  particles  derived  from  ethylene/air, 
surrogate/air,  m-xylene/air,  and  n-dodecane/air  flames  is  proposed  in  this  work  where: 

WMiglcnf  -sec)  =  1.2M02  %/>„„,  +Wa 

* T  (42) 


Wn  =kr 


*  T>n 


(43) 


ko2  is  the  intrinsic  reaction  rate  and  P 02  is  the  partial  pressure  in  atm.  ko2  can  be  expressed  as: 
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A*  T° 5  exp 


Ea 

~RT 


(44) 


where  A  is  the  pre-exponential  factor  ,  Ea  is  the  activation  energy,  R  is  the  gas  constant,  and  T  is 
the  reaction  temperature.  Equations  43  and  44  were  solved  numerically  in  Polymath  for 
conditions  where  oxidation  via  O2  dominated  close  to  the  burner  surface  (ethylene/air, 
Ooverall=0.8;  surrogate/air,  Ooverall  =  0.8,  1.14;  m-xylene/air,  ^overall  =  0.87,  1.14;  n- 
dodecane/air,  Ooverall  =  0.8,  1.14.  The  parameters  A,  n  and  Ea  that  best  fitted  the  experimental 
Ws  for  each  condition  are  presented  in  Table  5. 


Table  5.  Kinetic  parameters  that  produced  best  fitting  for  the  experimental  data. 


System 

A 

Ea,  kJ/mol 

n 

Ethy_0.8 

0.985 

161.50 

0.75 

Surrogate_0.8 

1.101 

140.37 

0.73 

Surrogate_1.14 

1.101 

126.61 

0.70 

m_Xylene_0.87 

1.101 

158.38 

0.77 

m_Xylene_1.14 

1.101 

105.43 

0.74 

n-Dodecane_0.8 

1.101 

126.96 

0.76 

n-Dodecane_1.14 

1.200 

100.40 

0.76 

Av  erase 

1.099 

131.38 

0.74 

The  averaged  values  presented  in  Table  5  were  used  to  express  the  overall  soot  oxidation  rate  as 
follows: 


Wtotal  (g!  cm2-  sec)  =  1 .29x1 02 


^ OH  p 

Vr  OH’ 


+  1  099*r0  5  exp 


UhAjUTmol) 

R*T 


(45) 
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This  expression  applies  for  equivalence  ratios  ranging  from  0.8  to  1.14,  O2  partial  pressures  from 
0.065  to  0.123  atm,  and  peak  temperatures  from  1450-1750  K  for  soot  oxidation  from 
ethylene/air,  surrogate/air,  m-xylene/air,  and  n-dodecane/air  flames.  Comparisons  of  predictions 
using  Equation  45  with  experimental  data  are  presented  in  Figures  45a  through  45g  for 
ethylene/air  (Figures  45a),  surrogate/air  (Figures  45b  and  45c),  m-xylene/air  (45d  and  45e)  and 
n-dodecane/air  (Figures  45f  and  45g)  premixed  flames.  Once  again,  the  effect  of  the  sampling 
probe  was  accounted  for.  It  is  also  interesting  to  note  that  the  activation  energy  is  approximately 
that  for  the  NSC  equation’s  A  sites  (see  Table  2)  which  is  30  kcal/mole  or  approximately  125 
kJ/mole.  The  rate  is  still  low  as  a  result  of  different  in  the  pre-exponential  factor,  although  the 
pre-exponential  is  a  function  of  surface  area.  As  previously  discussed,  an  increase  in  surface 
area  might  account  for  these  differences.  These  results  also  suggest  that  the  studies  were  in  the 
regime  of  temperatures  and  lower  reactivity  sites. 

Results  in  Figure  45  showed  a  good  agreements  between  experimental  soot  oxidation  rates 
determined  from  SMPS  measurements  and  predictions  using  the  proposed  Equation  45  which 
accounted  simultaneous  for  the  effect  of  temperature,  O2  and  OH*  concentration  on  the  soot 
oxidation  derived  from  different  sources  in  premixed  combustion  systems.  The  proposed 
equation  is  able  to  capture  the  major  features  of  the  experimental  results,  particularly  the  major 
contribution  of  O2  to  soot  oxidation  in  the  region  close  to  the  burner  surface  and  OH*  higher  in 
the  burner.  In  comparison  with  the  literature  (i.e.  Nienow,  et  al.  2005),  the  activation  energy  is 
within  the  range  of  others  as  is  the  order  in  oxygen.  Literature  values  tend  toward  100-150 
kcal/mol  with  an  order  of  approximately  0.7  in  oxygen. 

Changes  in  Soot  Morphology  and  Nanostructure  During  Soot  Oxidation 

Preliminary  TEM  and  HR-TEM  analysis  of  soot  samples  collected  thermophoretically  in  the  top 
burner  as  a  function  of  the  height  above  the  burner  for  ethylene/air  flame  under  a  lean  condition 
were  carried  out  to  compare  with  the  SMPS  data  and  obtain  insights  into  the  changes  on  soot 
morphology  and  nanostructure  during  soot  oxidation  and  soot  fragmentation.  Figure  46  shows 
the  results  for  the  ethylene/air  flame  (Ooverall  =  0.87).  Representative  soot  samples  at  different 
HAB  were  collected  using  a  more  heavily  sooting  flame  in  the  first  burner  (ethylene/air, <J>i  = 
2.8),  but  the  overall  equivalence  ratio  was  kept  constant  at  0.87.  PSDs  (see  Figures  34a  and  35a) 
for  this  condition  showed  the  presence  of  soot  fragmentation  close  to  the  burner  surface  and  soot 
burnout  mainly  higher  in  the  burner. 
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Figure  45.  Predictions  of  soot  oxidation  rates  using  Equation  24  for:  (a)  ethylene/air;  Ooverall  = 
0.8  (b)  and  (c)  surrogate/air,  Ooverall  =  0.8,  1.14;  (d)  and  (e)  m-xylene/air,  Ooverall  =  0.87, 
1.14  and  (f)  and  (g)  n-dodecane/air,  Ooverall  =  0.8,  1.14. 
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Figure  46.  TEM  pictures  for  ethylene/air  flame  with  ^overall  =  0.87.  (a)  At  HAB  =  2  mm  and 

(b)  At  HAB  =  4  mm. 


Figure  47  (left)  presents  the  results  for  the  evolution  of  the  oxidation  process  as  a  function  of  the 
HAB  (<0.5  mm,  2  and  4mm  for  this  condition).  As  seen  in  the  TEM  images,  close  to  the  burner 
surface,  most  of  the  particles  were  present  as  agglomerates.  As  HAB  increases  a  drop  in  the 
number  of  particles  is  shown  as  they  were  burned  out.  This  behavior  followed  the  trend  of  the 
SMPS  data  in  terms  of  particle  number  and  mass  concentration  decrease.  TEMs  of  collected 
samples  were  also  taken  during  the  oxidation  of  soot  derived  from  the  surrogate  mixture  in  air. 
TEM  images  (see  Figure  47,  right)  revealed  a  similar  behavior  as  compared  to  ethylene/air 
flames.  Agglomerates  of  particles  dominated  close  to  the  burner  surface.  As  the  oxidation 
process  continued,  particle  diameter  and  number  concentration  dropped.  These  results  were 
mainly  consistent  with  the  evolution  of  PSDs  as  a  function  of  the  HAB  in  terms  of  the  burnout  of 
particles. 

HR-TEM  picture  were  processed  and  analyzed  using  the  methodology  previously  described 
(Toth  et  al.  2010).  The  nanostructural  evolution  from  our  samples  collected  in  the  two-stage 
burner  along  with  other  samples  collected  from  fuel  rich  ethylene/air  and  benzene/air  (Echavarria 
et  al.  2009,  2011)  with  different  residence  times  and  the  extents  of  oxidation  is  shown  by 
monitoring  changes  in  distance  deviation  parameter,  O,  and  junction  parameter,  u.  Table  6  is  a 
summary  of  the  various  samples  and  their  names. 
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Figure  47.  TEM  pictures  for  ethylene/air  (left)  and  surrogate  (right)  at  HAB  =  0.5  mm,  2  mm, 
and  4  mm.  Note,  the  scale  of  the  top  image  is  different. 
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Table  6.  Summary  of  samples  collected 


Sample  name 

Fuel  type 

Burner  setup 

Heights  above  burner,  mm 

Benzene  flame 

Benzene/air 

Premixed,  fuel  rich 

5,  10,  15 

Ethylene  flame  1 

Ethylene/air 

Two  stage  premixed,  fuel  rich/fuel  lean 

0.5,  2, 4 

Ethylene  flame  2 

Ethylene/air 

Pie  mixed,  fuel  rich 

10,  15 

Surrogate  flame 

N-dodecane-m-xylene  mixture/air 

Two  stage  premixed,  fuel  rich/fuel  lean 

1,3,5 

The  parameters,  Q  and  v  were  calculated  for  four  fuels,  each  one  sampled  at  different  heights 
above  burner  (HAB).  As  HAB  increases  and  the  soot  gets  more  mature,  an  increase  in 
orderliness  can  be  expected  as  the  conventional  analysis  of  Alfe  et  al.  (2010)  demonstrated  for 
methane  soots.  A  few  examples  of  the  analyzed  images  can  be  seen  in  Figure  48.  It  can  be 
generally  stated,  that  among  the  soots  collected,  ethylene  flames  produced  the  most  unordered 
structures,  while  surrogate  flames  produced  the  most  ordered  structures,  exhibiting  high 
concentric  symmetry  at  higher  HABs.  The  HR-TEM  micrographs  were  analyzed  by  both 
algorithms.  In  the  case  of  the  determination  of  f2,  the  Euclidean  distance  set  of  each  image  was 
extracted,  and  added  to  datasets  representing  each  sample.  These  datasets  were  then  evaluated,  to 
determine  Q.  The  average  number  of  distances  extracted  from  one  512  x  512  image  is  ~  50,000, 
and  several  micrographs  were  analyzed  for  each  sample,  v  was  calculated  for  each  micrograph 
and  average  values  were  produced  for  different  samples. 

First,  a  conventional  analysis  of  these  images  was  conducted  for  the  sake  of  comparison,  using  a 
methodology  based  on  the  technique  published  by  Palotas  et  al.  (1996).  The  measured 
parameters  of  the  fringe  images  were  the  fringe  lengths,  [nm]  -  measured  as  proposed  by 
VanderWal  et  al.  (2010),  i.e.  the  length  of  the  skeleton  of  the  fringe  -  and  interlayer  distances, 
[nm],  measured  as  proposed  by  Palotas  et  al.  (1996).  The  results  of  this  statistical  analysis  are 
presented  in  Figure  49.  As  shown  in  Figure  49,  conventional  statistical  analysis  produces  results 
based  on  which  the  reliable  differentiation  of  similar  structures  is  difficult.  The  values  of  the 
interlayer  distances  show  insignificant  differences  between  samples  and  there  is  no  apparent 
trend  in  the  evolution  of  fringe  lengths.  Apart  from  this,  the  inherently  high  standard  deviation  of 
the  results  provided  by  this  low-level  image  analysis  algorithm  makes  the  differentiation 
impossible,  based  on  mean  values  alone. 

The  results  for  Q  and  v  are  presented  in  Figure  50.  The  trends  of  the  two  parameters  are  quite 
obvious  as  both  decreased  with  increasing  HABs  for  all  four  samples,  meaning  that  an  increasing 
magnitude  of  order  could  be  detected  in  the  soot  structures.  Interestingly  enough,  an  even 
stronger  correlation  was  found  between  Q  and  v.  The  highest  degree  of  structural  evolution 
happened  in  the  soots  from  the  ethylene  flames,  especially  at  low  HABs.  The  rather  small  change 
in  Q  measured  in  surrogate  soots  can  be  explained  by  the  fact  that  these  soots  were  found  to  be 
fairly  well  ordered  even  at  low  HABs,  however  a  huge  drop  in  v  indicates  some  extent  of 
simplification  in  the  texture  (notice  the  highly  developed  concentric  rings  in  Figure  48  D). 
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Figure  48.  Some  examples  of  the  analyzed  images.  A,  Ethylene  flame  1,  HAB=0.5  mm.  B, 
Ethylene  flame  1,  HAB=4  mm.  C,  Surrogate  flame,  HAB=1  mm.  D,  Surrogate  flame,  HAB=5 
mm.  E,  Benzene  flame,  HAB=5  mm.  F,  Benzene  flame,  HAB=15  mm. 
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Height  above  burner,  mm 

Figure  49.  Results  of  a  conventional  statistical  image  analysis  method  applied  to  our  samples. 
Some  HAB  values  were  shifted  to  improve  visibility.  For  exact  values  see  Table  6. 
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Figure  50.  The  measured  values  of  Q  and  o  for  the  fuel  samples.  Lines  are  guides 

to  the  eye. 


Feasibility  of  Soot  Fragmentation  due  to  Internal  Burning 
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The  process  of  soot  oxidation,  as  evaluated  from  detailed  and  spatially  resolved  measurements  of 
soot  size  distribution,  flame  temperature,  gas-phase  composition  (O2  and  OH*),  and  soot 
morphology  and  nanostructure  has  provided  critical  information  on  the  kinetics  of  soot  oxidation 
via  O2  and  OH*.  A  phenomenon  where  particle  concentration  in  the  ultrafine  (Dp  <  lOnm)  mode 
(see  Figures  3 1  and  40)  increases  with  HAB  during  oxidation  has  been  associated  to  the  break  up 
or  fragmentation  of  the  particles  which  exist  as  agglomerates  rather  than  individual  particles  (see 
Figure  47).  A  major  difference  between  the  results  in  this  work  and  that  of  Neoh  et  al.  (1981) 
was  that  the  observable  increase  of  nanoparticles  in  the  ultrafine  mode  occurred  mainly  at  low 
soot  burnouts  for  this  work  and  at  high  conversions  for  Neoh  et  al.  At  high  burnout  percentages, 
this  behavior  is  believed  to  be  consistent  with  percolation  theory  which  predicts  fragmentation 
with  varying  pore  structure  (Kerstein  et  al.,  1985).  For  the  present  work,  it  appears  as  though  the 
fragmentation  involved  the  bridges  between  primary  particles  and  occurred  at  low  overall  soot 
burnout  (<  15%). 

Figure  51  shows  a  schematic  representation  of  the  fragmentation  process  at  low  and  high  burnout 
percentages,  particularly  under  overall  lean  conditions  where  O2  can  also  contribute  to  soot 
oxidation  higher  in  the  burner.  According  to  Neoh  and  coworkers,  fragmentation  can  be 
observed  at  high  burnout  percentages  (see  Figure  51,  right  side)  when  high  internal  surface  area 
or  high  porosity  has  been  developed  in  the  particle  surface.  Neoh  et  al.  (1981)  attributed  this 
phenomena  to  O2  penetrating  the  soot  nanostructure,  which,  in  turn,  caused  internal  burning  and 
loss  of  the  connectivity  between  the  carbonaceous  phases  inside  the  particle.  This  situation 
resulted  in  the  fragmentation  or  break  up  of  primary  particles.  We  have  attributed  the 
fragmentation  at  low  burnout  to  both  the  breakup  of  bridges  (or  necks)  between  primary 
particles  and  the  removal  of  small  fragments  from  the  edges  of  the  primary  particles  surface  as 
illustrated  in  Figure  51,  left  side. 
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Figure  51.  Hypothesis  regarding  the  role  of  fragmentation  during  soot  oxidation  at  low  soot  burnout  (left)  and  higher  soot  burnout 


Conditions  that  favor  this  mechanism  can  be  identified  from  calculations  of  the  effectiveness 
factor  (rj)  for  O2  and  OH*.  The  effectiveness  factor  is  the  ratio  between  the  actual  oxidation  rate 
and  the  oxidation  rate  assuming  that  the  whole  internal  surface  was  exposed  to  the  reactant  and 
that  the  concentration  of  the  oxidizer  species  and  temperature  are  uniform  throughout  the 
particle,  an  idea  of  the  importance  of  internal  burning.  If  the  effectiveness  factor  is  close  to  1 ,  it 
is  an  indication  of  penetration  into  the  soot  nanostructure.  In  contrast,  if  the  effectiveness  factor 
is  close  to  zero  it  is  an  indication  of  little  gas  penetration  (Neoh,  1981). 


The  effectiveness  factor  can  be  expressed  as  (Neoh,  1981): 
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The  reaction  constant  ksi  was  calculated  using  both  Equation  45,  which  calculated  the  rate  of 
oxidation,  and  assuming  a  first  order  reaction  rate  (Neoh,  1981).  Results  were  similar  for  both 
conditions.  In  this  equation,  sv  is  the  pore  surface  (cm  /cm  )  and  equivalent  to  the  porosity 
divided  by  the  pore  radius.  The  effective  diffustion  coefficient,  Deff  (cm2/s),  is  calculated  from: 

Deff  =  0^).(1/Dg  +1  !Dkyx  (48) 

where  the  tortuosity  factor,  x,  was  assumed  to  be  2.  The  Knudsen  diffuction  coefficient,  Dk  is: 


D ^  =  9700/^  4t  1  Mwcm  ^  /  s 


(49) 


Values  for  the  tortuosity  and  pore  radius  represented  average  values  typically  found  in  the 
literature  for  different  types  of  soot  and  also  extracted  from  our  transmission  electron 
microscopy  images.  A  pore  radius  in  the  range  from  5  A  to  20  A  can  be  taken  to  represent  actual 
pores  in  individual  particles  and  the  spacing  between  individual  particles  in  the  agglomerates. 

Tables  7  and  8  show  the  results  of  r|  obtained  at  temperature  and  mean  particle  size  (averaged  in 
the  region  close  to  the  burner  surface)  conditions  representing  the  region  below  the  tube  bundle 
and  close  to  the  top  burner  surface  before  fragmentation.  The  objective  was  to  evaluate  the  extent 
of  O2  and/or  OH*  penetration  into  the  soot  porous  network  which  was  expected  to  contribute  to 
internal  burning.  For  the  ethylene/air  flames  with  Ooverall  =  0.8  (Table  7),  r)o2  was  1  or  close  to 
1  relative  to  the  smaller  values  obtained  for  r|0H*-  This  was  an  indication  that  the  O2  molecule 
would  be  able  to  penetrate  into  the  porous  network  to  a  higher  extent  than  OH*.  Therefore, 
internal  burning  might  occur  mainly  by  the  action  of  O2.  Similar  results  were  obtained  for  the 
surrogate  mixture  (see  Table  8).  At  higher  burnout  percentages,  the  degree  of  penetration  by 
OH*  becomes  significant  (Neoh,  1981);  however,  internal  burning  due  to  OH*  is  not  expected  to 
have  a  major  impact  due  to  its  higher  intrinsic  reactivity,  and  oxidation  via  OH*  would  be 
expected  to  take  place  mainly  in  the  particle  surface. 
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Table  7.  Effectiveness  factor  calculation  for  the  ethylene/air  system  with  ^overall  =  0.8 


Effectiveness  factor(h),  02 


373  K 

600  K 

1403  K 

1510,  K 

re  =  5  A 

1.00000 

1.00000 

1.00000 

1.00000 

re  =  10  A 

0.99998 

0.99998 

0.99995 

0.99991 

re  =  15  A 

0.99999 

0.99999 

0.99998 

0.99996 

re  =  25  A 

1.00000 

1.00000 

0.99999 

0.99998 

Effectiveness  factor(h),  OH* 

373  K 

600  K 

1403  K 

1510, K 

re  =  5  A 

0.1056 

0.1056 

0.1056 

0.1056 

re  =  10  A 

0.2032 

0.2032 

0.2032 

0.2032 

re  =  1 5  A 

0.2927 

0.2927 

0.2927 

0.2927 

re  =  20  A 

0.3743 

0.3743 

0.3743 

0.3743 

Mean  Particle  Diameter  (0-2.5  mm)  = 

43  nm 

Table  8.  Effectiveness  factor  calculation  for  the  surrogate  mixture/air  system  with  ^overall  =  0.8 


Effectiveness  factor(h),  02 


373  K 

600  K 

927K 

1500K 

re  =  5  A 

1.000 

1.000 

1.000 

1.000 

re  =  10  A 

0.999 

0.998 

0.998 

0.997 

re  =  15  A 

0.999 

0.999 

0.999 

0.999 

re  =  25  A 

1.000 

1.000 

0.999 

0.999 

Effectiveness  factor(h),  OH* 


373  K 

600  K 

927K 

1500K 

re  =  5  A 

0.038 

0.038 

0.038 

0.038 

re  =  10  A 

0.076 

0.076 

0.076 

0.076 

re  =  1 5  A 

0.112 

0.112 

0.112 

0.112 

re  =  20  A 

0.1476 

0.1476 

0.1476 

0.1476 

Mean  Particle  Diameter  (0-2.5  mm)  =  120  nm 


Another  indication  of  internal  burning  was  extracted  from  experimental  measurements  of  the 
specific  soot  surface  area  of  samples  collected  isokinetically  using  the  methodology  described  in 
the  materials  and  methods  section.  As  mention  previously,  soot  surface  area  was  determined 
using  CO2  adsorption  at  298  K  in  a  thermogravimetric  analyzer  (TGA).  Samples  collected  at  the 
second  burner  surface  (HAB  =  0  mm)  and  at  2.5  mm  above  second  burner  surface  showed  a 
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small  increase  in  soot  surface  area  for  both,  soot  derived  from  ethylene/air  (Ooverall  =  0.8)  and 
from  the  surrogate/air  (^overall  =  0.8)  flames.  For  ethylene  air,  soot  surface  area  increased  from 
282  ±  35  (0mm)  to  494  ±  27  (2.5  mm)  m2/gsoot.  For  the  surrogate/air  system,  the  surface  area 
increased  from  ~196  ±  18  to  ~285  ±  31  m2/gsoot  in  the  region  between  0  to  2.5  mm  respectively. 
Even  though  these  changes  are  relatively  small  in  both  flames,  the  increase  in  surface  area  was 
an  indication  of  internal  burning  that  would  be  expected  to  favor  the  removal  of  fragments  from 
edges  of  primary  particles  as  illustrated  in  Figure  5 1  (left  side). 

TEM  and  HR-TEM  analysis  provided  evidence  of  the  presence  of  the  type  of  structures  that 
might  favor  the  mechanisms  described  above.  For  example,  the  presence  of  chain  like 
agglomerates,  bridges  or  necks  cementing  primary  particles,  and  the  spherical  characteristics  of 
the  primary  particles  close  to  the  second  burner  surface  (0-2.5  mm).  Higher  in  the  burner  (HAB 
=  5  mm),  the  presence  of  fragments  without  a  well  defined  form  was  also  observed,  and  their 
presence  was  associated  to  the  breakup  of  primary  particles.  Additional  HR-TEM  analysis 
allowed  establishing  a  correlation  between  the  change  in  the  diameter  of  the  necks  between 
primary  particles  and  HAB.  Although  a  more  deep  statistical  analysis  of  the  HR-TEM  picture 
was  difficult  to  perform  higher  in  the  burner  because  of  limited  sample,  Figure  52  shows  that  the 
diameter  of  the  necks  decreased  with  increasing  HAB,  and  it  would  be  expected  to  favor  the 
breakup  of  the  necks  cementing  primary  particles.  This  mechanism  and  the  removal  of  fragments 
from  the  edges  of  primary  particles  in  the  region  close  to  the  burner  surface  seem  to  be  an  initial 
support  for  the  observable  increase  of  nanoparticles  in  the  ultrafine  mode  from  SMPS 
measurements. 
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Figure  52.  Change  of  the  diameter  in  the  necks  or  bridges  cementing  primary  particles  with 
HAB.  Analysis  was  performed  for  20  pictures  for  each  height. 


CONCLUSIONS  AND  IMPLICATIONS 

This  project  focused  on  the  mechanisms  of  soot  oxidation  by  oxygen  for  jet  fuels  used  in  military 
engines.  The  reactions  with  oxygen  are  particularly  important  for  gas  turbines,  which  operate  at 
high  excess-air  rates.  The  examination  of  the  role  of  fuel  components  and  internal  surface  on 
soot  oxidation  rates  was  emphasized  since  the  allowance  for  these  effects  on  the  oxidation  rate 
are  not  currently  well  characterized.  Examination  of  this  additional  aspect  of  oxidation  kinetics 
enhanced  fundamental  combustion  science  and  over  the  long  term  will  enable  military  gas- 
turbine  engines  to  continue  to  meet  their  performance  and  operating  requirements  with  reduced 
PM  emissions. 

Soot  oxidation  was  studied  by  following  the  evolution  of  particle  size,  particle  number 
concentration,  flame  temperature,  gas-phase  composition,  particle  surface  area,  and  particle 
nanostructure  and  morphology  during  the  oxidation  of  soot  derived  from  ethylene/air, 
surrogate/air,  m-xylene/air,  and  n-dodecane/air  premixed  flames  in  a  two-stage  burner  where 
soot  formation  mechanisms  were  isolated  from  the  oxidation  steps.  Soot  was  burned  under  fuel- 
lean  and  slightly  fuel-rich  conditions.  In  ethylene/air  flames  ,  the  analysis  of  the  results  in  terms 
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of  PSDs  and  number  concentration  measurements  showed  the  possibility  of  break  up  or 
fragmentation  of  the  bigger  particles,  as  evidenced  by  a  significant  increase  of  particles  in  the 
ultrafine  mode  (Dp<10nm)  and  decrease  in  particle  mean  diameter.  The  extent  of  the 
fragmentation  decreased  with  increasing  the  overall  equivalence  ratio  from  0.8  to  0.94,  and  it  did 
not  take  place  for  rich  flames  (Ooverall  =  1.07  and  1.14).  Analysis  of  the  results  for  ethylene/air 
flames  in  terms  of  the  evolution  of  major  gas-phase  species,  particularly  soot  oxidizer  species  in 
flames  like  O2  and  OH*  provided  information  about  the  contributions  of  O2  and  OH*  on  soot 
fragmentation  and  soot  burnout.  For  lean  conditions  (Ooverall  =  0.8  and  0.87),  O2  oxidation  was 
predominant  close  to  burner  surface,  and  appeared  to  favor  the  fragmentation  of  particles. 
Downstream,  OH*  had  a  major  impact,  and  soot  oxidation  rates  increased.  For  the  rich  flames 
(Ooverall  1.07  and  1.14),  soot  oxidation  rates  were  higher  close  to  burner  surface  which  was 
attributed  to  the  faster  consumption  of  O2  and  earlier  formation  of  OH*. 

Results  during  the  oxidation  of  soot  particles  derived  from  the  JP-8  surrogate  mixture  and  its 
pure  components  (m-xylene  and  n-dodecane)  showed  a  similar  behavior  under  lean  conditions 
(Ooverall  =  0.8,  0.87)  relative  to  ethylene/air  flames.  However  the  fragmentation  phenomenon 
was  also  observed  to  occur  under  slightly  rich  conditions  (Ooverall  =  1.14),  and  the  extent  of 
fragmentation  was  as  high  as  or  even  higher  than  the  one  observed  under  lean  conditions.  This 
behavior  was  correlated  with  the  evolution  of  gas-phase  species  that  showed  significant  changes 
in  the  evolution  of  soot  oxidizer  species  after  a  prolonged  induction  period  where  high 
concentrations  of  O2  and  a  high  temperature  environment  favored  the  fragmentation  process. 

A  major  difference  between  the  results  in  this  work  and  those  of  Neoh  et  al.  (1981)  was  that  the 
fragmentation  occurred  mainly  at  low  soot  burnouts  for  this  study  and  at  high  conversions  for 
Neoh  et  al.  (1981).  For  the  present  work,  it  appears  as  though  the  fragmentation  involved  the 
bridges  between  primary  particles  and  occurred  at  low  overall  soot  burnout. 

The  fragmentation  process  was  associated  to  O2  penetration  into  the  soot  porous  network,  which 
may  have  caused  internal  burning  and  resulted  in  the  breakup  of  both  the  bridges  (necks) 
between  primary  particles  and  the  rupture  of  primary  particles  that  might  occur  close  to  the  top 
burner  surface  and/or  higher  elevations.  The  feasibility  of  these  types  of  mechanisms  was 
supported  by  calculations  of  the  effectiveness  factor  and  experimental  measurements  of  soot 
surface  area  and  morphology.  The  effectiveness  factor  calculations  in  ethylene/air  and 
surrogate/air  flames  suggested  that  the  O2  molecule  would  be  able  to  effectively  penetrate  the 
soot  nanostructure  in  the  region  close  to  the  burner  surface.  Internal  burning  was  confirmed  by 
the  increase  in  soot  surface  area  in  that  region  (0  to  2.5  mm).  TEM  and  HR-TEM  analysis 
revealed  the  morphological  and  nanostructural  characteristics  of  soot  particles  that  can  support 
the  breakup  of  the  bridges  cementing  primary  particles  and  the  removal  of  fragments  from  the 
particles  edges  due  to  internal  burning. 

TEM  and  HR-TEM  micrographs  also  showed  insights  into  the  changes  in  soot  morphology  and 
nanostructure  that  occur  during  soot  oxidation  and  fragmentation.  TEM  pictures  supported  the 
burnout  of  particles,  while  analysis  of  HR-TEM  pictures  showed  difference  of  soot  derived  from 
ethylene/air  and  the  surrogate/air  flames.  A  sensitive  image  analysis  technique  has  been 
developed  for  the  representative  quantification  of  soot  nanostructure  orderliness  utilizing  image 
processing  techniques.  The  technique  is  robust,  since  the  evaluated  number  of  elements  in  the 
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measured  datasets  is  large,  and  produces  consistently  reproducible  parameters  with  a  resolution 
high  enough  to  show  differences  in  soot  structure.  The  method  has  been  applied  to  soot  samples 
of  various  fuels,  collected  from  different  burner  setups  at  different  HAB  in  order  to  monitor  their 
nanostructural  evolution.  In  all  cases,  decreasing  tendencies  were  found  in  the  values  of  the 
distance  deviation  and  junction  parameters  as  the  HAB  increased  suggesting  increasing 
orderliness  in  soot  structure  with  its  age.  Based  on  the  observation  of  Q  ,  the  highest  degree  of 
structural  evolution  happened  in  the  soots  from  the  ethylene  flames,  especially  at  low  HABs.  The 
rather  small  change  in  Q  measured  in  surrogate  soots  can  be  explained  by  the  fact  that  these 
soots  were  found  to  be  fairly  well  ordered  even  at  low  HABs.  However,  a  huge  drop  in  u 
indicates  some  extent  of  simplification  in  the  texture  as  reflected  in  the  highly  developed 
concentric  rings  in  soot  derived  from  the  surrogate  flame. 

The  experimental  data  were  compared  to  typical  models  used  in  the  literature  and  they  produced 
a  good  agreement  with  the  experimental  results.  A  kinetic  expression  that  accounted  for  changes 
in  the  temperature,  and  O2,  and  OH*  concentration  during  the  oxidation  of  soot  produced  from 
different  sources  was  proposed.  The  proposed  equation  was  able  to  reproduce  the  major  features 
of  the  experimental  results  and  cover  wide  range  of  conditions  typical  of  practical  combustion 
systems.  In  addition  the  activation  energy  and  oxygen  partial  pressure  order  agreed  within  the 
range  presented  in  the  literature.  The  pre-exponential  was  different,  but  this  is  likely  a  result  of 
changes  in  the  surface  area.  Even  though  the  nanostructure  changed  with  HAB  for  the  ethylene 
and  surrogate,  it  was  not  apparent  that  this  change  in  nanostructure  affected  the  rate  of  oxidation 
or  that  our  data  accounted  for  such  changes. 

The  design  of  the  experimental  and  modeling  methodologies  on  soot  oxidation  in  premixed 
flames  allowed  for  identifying  the  conditions  where  soot  fragmentation  and  soot  burnout  can  be 
favored.  New  insights  into  the  process  of  soot  fragmentation  were  recognized  by  showing  that 
fragmentation,  depending  upon  the  experimental  conditions,  can  occur  at  lower  soot  burnout 
percentages  and  under  fuel-rich  conditions.  These  results  are  expected  to  contribute  to  the  current 
knowledge  on  the  oxidation  of  soot  particles.  This  study  also  opens  new  paths  to  evaluate 
different  strategies  to  reduce  the  emissions  of  soot  into  the  atmosphere. 
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